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IV 


ABSTRACT 


Mathematicians  have  used  quaternions  for  about  a  hundred  years.  Today  they  are 
an  important  part  of  computer  graphics  and  simulation  systems.  This  thesis  takes  an 
analytical  approach  to  quaternions  by  using  them  in  the  construction  of  a  virtual  human 
for  sourceless  Magnetic  Accelerometer  Rate  Sensor  (MARG)  body  tracking  system. 

Virtual  citizens  will  be  a  reflection  of  our  personalities  in  cyberspace.  Prophecies 
say  they  may  take  control  in  the  virtual  world  and  govern  themselves  too.  One  of  the 
objectives  of  this  thesis  is  to  design  a  seamless  and  realistic  humanoid  from  laser  scan 
data  clouds.  This  humanoid  will  be  compatible  with  motion  capture  systems  and 
networked  virtual  environments. 

Second  objective  of  this  thesis  is  to  search  for  the  answers  related  to  the  optimal 
real-time  representation  of  an  articulated  virtual  human,  maintaining  a  high  level  of 
visual  fidelity  within  networked  cyberspace.  While  visual  detail  and  fidelity  have  been 
and  will  continue  to  be  a  major  ongoing  interest  within  the  computer  graphics 
community,  the  idea  of  sourceless  body  tracking  is  still  in  its  early  stages.  MARG  body 
tracking  is  one  of  the  successful  approaches  to  body  tracking  systems.  This  thesis 
proposes  a  networked  quaternion  based  real-time  virtual  human  interface  for  the  MARG 
body  tracking  system. 
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I.  INTRODUCTION 


A.  MOTIVATION 

Since  the  earliest  day  of  the  art  form,  animators  have  observed  the  movement  of 
real  ereatures,  in  order  to  ereate  animated  motion.  Sometimes  this  simply  takes  the  form 
of  an  artist  carefully  observing  nature  for  inspiration.  Another  proeess  is  to  transfer  the 
movement  from  reeording  of  the  movement  to  the  animated  objeets.  The  earliest 
meehanism  for  doing  this  was  the  Rotoseope,  a  device  that  projeeted  frames  of  film  onto 
the  animator’s  workspace,  providing  the  animator  with  a  guide  for  their  drawings. 

Today  eomputer  power  and  improvements  in  miero  machine  teehnology  allow 
researchers  to  design  special  sensors  in  order  to  traek  the  human  body.  The  magnetie, 
Aeeelerometer,  Rate  sensor  (MARG)  based  body  tracking  system  is  an  ongoing  projeet  at 
the  Naval  Postgraduate  Sehool  and  is  a  good  example  of  this  technology.  The  MARG 
body  traeking  system  does  not  depend  on  an  external  souree  for  body  traeking  and  this 
provides  a  remote  body  traeking  eapability  to  the  system  with  a  speeial  position-traeking 
deviee.  This  thesis  eoneentrates  on  using  outputs  of  a  MARG  system  for  a  virtual  human. 


B,  PROBLEM  STATEMENT 

Even  though  there  is  a  big  inerease  in  eentral  proeessor  and  graphics  processor 
units,  it  is  still  too  diffieult  to  animate  full  human  body  motion  in  real  time.  There  had 
been  a  previous  work  to  aehieve  this  goal  by  using  segmented  human  body 
[DUTTON2001].  The  segmented  human  body  has  a  joint  break  problem  when  it  moves. 
This  thesis  foeused  on  this  problem  and  tried  to  answer  the  question  of  how  to  ereate  a 
seamless  skin  model  with  segmented  skeleton,  whieh  uses  quaternion  information  from 
the  MARG  body  tracking  system. 

C.  PROPOSED  SOLUTION 

As  stated  in  the  problem  statement  this  thesis  looks  for  a  method  to  eonstruet  a 
virtual  human.  [DUTTON2001]  proposed  a  solution  to  design  a  virtual  human  from  laser 


1 


scans  by  using  H-Anim  1.1  specification.  This  research  follows  the  same  way  with  a 
different  approaeh.  Instead  of  segmented  skeleton  structure,  this  thesis  used  a  skin  model 
with  a  skeleton  structure,  which  is  an  application  of  H-Anim  2000  standard.  The  virtual 
human  will  have  the  ability  to  position  itself  by  using  Global  Position  System  data  and  to 
orient  its  segments  by  using  MARG  sensor  quaternion  data 


D,  DESIGN  CONSIDERATIONS  AND  CONSTRAINTS 

There  are  three  constraints  for  this  research.  The  first  constraint  is  related  to  the 
real-time  rendering  capabilities  of  computer  graphics  systems.  As  the  level  of  detail  for 
laser  scans  increases  ,the  speed  of  real-time  rendering  decreases.  So  this  trade  off 
between  realistic  appearance  and  real-time  rendering  has  leaded  this  research  to  search 
for  an  optimum  solution. 

The  second  constraint  is  the  number  and  current  capabilities  of  MARG  sensors. 
There  are  only  three  MARG  sensors  finished  right  now.  This  prevents  experimentation 
with  full  body  articulation.  A  simulated  quaternion  data  broadcasting  system  was 
designed  because  of  the  lack  of  real-time  networking  for  processed  data. 

Finally  it  is  noted  that  currently  there  is  no  eomputer  graphies  system  that  uses 
quaternions  for  graphical  transformation  purposes.  A  quaternion  to  axis  angle 
transformation  class  is  written  to  simulate  the  concept. 


E,  GOALS 

The  objective  of  this  thesis  is  to  analyze  the  quaternion  rotations  in  detail  and  to 
construct  a  virtual  human  interface  for  a  quaternion  based  MARG  motion  eapture  system. 
Designing  a  realistie  and  real-time  virtual  human  interface  has  three  steps.  First  step  is 
the  construction  of  a  realistic  appearance.  This  is  achieved  by  using  laser  sean  data  eloud. 
The  seeond  step  is  segmenting  this  static  virtual  human  data  cloud  for  animations  by 
using  H-Anim  2.0  specification.  The  final  step  is  interfaeing  virtual  human  and  MARG 
body  tracking  system  through  the  networked  virtual  environment. 
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F.  ORGANIZATION 

Chapter  II  of  this  thesis  provides  some  fundamental  information  about 
quaternions.  Espeeially  for  definition  of  quaternions,  Sir  William  Hamilton’s  method, 
which  starts  with  vector  algebra,  is  used.  Also  in  this  chapter  some  basic  information 
about  MARG  Body  Tracking  System  and  Virtual  Humans  is  given.  Chapter  III  provides 
an  overview  of  modeling  a  virtual  Human  from  Laser  Scans  and  its  deficiencies.  Chapter 
IV  discusses  the  application  of  quaternion  based  MARG  sensors  in  Virtual  Environments. 
Chapter  V  is  the  summary  and  the  conclusions  from  the  experience.  Because  the  thesis  is 
part  of  an  ongoing  project ,  this  chapter  also  proposes  some  future  topics  for  research. 
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II.  BACKGROUND 


A.  QUATERNIONS 


1.  Elements  Of  Quaternion 

The  quotient  of  two  vectors  is  called  a  Quaternion,  .[HAMILTON1899] 

Let  a  and  p’  (Fig.l)  be  two  veetors  drawn  from  O  and  O’  respectively  and  not 
lying  on  the  same  plane;  and  let  their  quotient  be  designated  in  the  usual  way  by  a  /  P’  . 


Figure  1.  Quaternions  as  quotient  of  two  vectors 

Whatever  their  relative  positions,  we  may  always  conceive  that  one  of  these 
vectors,  as  P’,  may  be  moved  parallel  to  itself  so  that  the  point  O’  shall  move  over  the 
line  00’  to  O.  The  vectors  will  then  lie  in  the  same  plane  since  neither  the  direction  nor 
the  length  of  P’  has  been  changed  during  this  parallel  motion,  we  have  P  =  P‘.  And  the 
quotient  of  any  two  vectors,  a  and  P’  will  be  the  same  as  two  equal  co-initial  vectors  as  a 
and  p.  We  are  than  to  determine  the  ratio  a/p  in  which  a  and  P  lie  in  the  same  plane  and 
have  a  common  origin  O. 
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Whatever  the  nature  of  this  quotient,  we  are  to  regard  it  as  some  faetor  whieh 
operating  on  the  divisor  produees  the  dividend  i.e.  eauses  p  to  eoineide  with  a  in 
direetion  and  length,  so  that  if  this  quotient  is  q,  we  shall  have  by  definition. 


q.p  =  a 


(Eq.l) 


when 


a 


(Eq.  2) 


If  at  the  point  O'  we  suppose  a  veetor  O'C  =  7  to  be  drawn,  not  parallel  to  the 
plane  AOB,  and  that  this  veetor  be  moved  as  before,  so  that  O’  falls  at  O,  the  plane  whieh 
after  this  motion,Y  will  eoineide  on  a,  will  differ  from  the  plane  of  a  and  P,  so  that  if  the 
quotient;  a/p  =  q'  ,  q  and  q'  will  differ  beeause  their  planes  differs.  Henee  we  eonelude 
that  the  quotients  q  and  q'  eannot  be  the  same  if  a  ,  P  and  y  are  not  parallel  to  one  plane, 
therefore  that  the  position  of  the  plane  of  a  and  P  must  enter  into  our  eoneeption  of 
quotient  q. 

Again,  if  y  is  a  veetor  O'C,  parallel  to  the  plane  AOB,  but  differing  as  a  veetor 
from  P’,  then  when  moved,  as  before,  into  the  plane  AOB,  it  will  make  with  a  an  angle 
other  than  BAO.  Henee  the  angle  between  a  and  P  must  also  enter  into  our  eoneeption  of 
q.  The  direetion  of  the  angle  is  also,  ealled  direetion  of  rotation  another  important 
eontribution  to  q. 

Einally,  if  the  lengths  of  P  and  y  differ,  then  a/p  =  q  will  still  differ  from  a/y  =  q'. 
Therefore  the  lengths  of  veetors  must  also  enter  into  eoneeption  of  q. 

We  have  found  the  quotient  q,  regarded  as  an  operator,  whieh  ehanges  P  into  a, 
depending  upon  the  plane  of  veetors,  the  angle  between  them  and  the  ratio  of  their 
lengths.  Sinee  two  angles  are  requisite  to  fix  a  plane,  it  is  evident  that  q  depends  on  four 
elements  and  performs  two  distinet  operations.  Eirst  operation  is  named  as  stretehing  of 
P,  as  to  make  it  the  same  length  as  a  and  the  seeond  operation  is  turning  of  p,  as  to  eause 
it  to  eoineide  with  a  in  direetion. 
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Of  the  four  elements  of  quaternion; 

•  Turning  operation  depends  on  three  angles 

Two  angles  to  fix  the  plane  of  rotation 
One  angle  to  fix  the  amount  of  rotation 

•  Stretehing  depends  on  the  ratio  of  the  veetor  lengths 


A  quaternion  is  a  eomplex  quantity,  whieh  is  deeomposable  into  two  faetors.  One 
of  them  stretehes  and  shortens  the  veetor  divisor  in  order  to  provide  the  same  length  of 
veetor  dividend.  This  is  ealled  Tensor  of  the  quaternion.  Other  part  serves  to  turn  the 
divisor  for  ooineidenee  into  dividend  veetor,  whieh  is  ealled  as  Versor  of  the  quaternion 
by  Sir  William  Hamilton. [HAMILTONl 899]  These  factors  are  denoted  as  T(q)  and  U(q) 
respectively.  So  q  maybe  written 

q=  T(q)  .  U(q).  (Eq.3) 

2,  Properties  of  a  Quaternion 

a.  Equality  of  Quaternions 

By  using  previous  considerations  we  can  define  the  equality  of 
quaternions.  If  the  plane  of  a  and  P  is  moved  parallel  to  itself  or  if  the  angle  AOB  (Fig 
1)  ,  remaining  constant  in  magnitude  and  measured  in  the  same  direction,  is  rotated  about 
an  axis  through  o  perpendicular  to  the  plane;  or  the  absolute  lengths  of  a  and  P  vary  so 
that  their  ratio  remains  constant,  q  will  remain  same. 

Simply,  define 

a  /  P  =  q  and  a'  /  P'  =  q' 

q  =  q’  if  and  only  if 

•  The  vector  lengths  are  in  the  same  ratio  and 

•  The  vectors  are  in  the  same  or  parallel  planes,  and 
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•  The  vectors  make  with  each  other  the  same  angle  both  in  magnitude  and 
direction. 

We  should  recall  that  again  that  q  is  treated  as  an  operator  that  converts  P 
into  a  by  qP  =  a. 

b.  Quadrantal  Versors 

As  quaternion  q  operates  on  vector  P  it  involves  not  only  as  rotation  angle 
but  also  as  the  direction  of  this  rotation  angle.  This  leads  us  to  the  concept  of  a  reference 
frame.  By  positive  direction  we  assume  the  simple  right  hand  rule. 


Figure  2.  Quadrantal  Versors 


Let  i,  j  ,k  (Fig.  2)  represent  unit  vectors  at  right  angles  to  each  other.  The 
effect  of  any  unit  vector  acting  as  a  multiplier  upon  another  at  right  angles  to  it  has  been 
defined  to  be  the  turning  of  the  latter  in  a  positive  direction  in  a  plane  perpendicular  to 
the  operator  or  multiplier  through  an  angle  of  90°.  Thus  i,  operating  on  j,  produce  k.  This 
operation  is  called  multiplication  and  is  expressed  usually; 

ij  =  k.  (Eq.4; 

Another  approach  is  a  result  of  quotient  representation; 


k/j  =  i. 
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(Eq.5) 


It  has  to  be  pointed  out  that  this  multiplication  is  not  that  of  algebra;  it  is  a 
revolution.  Neither  k  in  Equation  4  is  a  numerical  product  nor  i  in  Equation  5  is  a 
numerical  quotient.  This  kind  of  multiplication  is  called  geometric. 

The  following  table  may  be  obtained  in  accordance  with  above. 


ij  =  k 

k/j  =  i 

II 

i/k  =j 

ki  =j 

J/i  =  k 

Ji  =  -k 

-k/i  =J 

kj  =  -I 

-i/j  =  k 

II 

-j/k  =i 

Table  1.  Products  of  Quadrinomial  Versors 

Because  the  effect  of  i,  j,  k,  is  a  rotation  from  one  direction  to  another  with 
an  angle  of  90°  they  are  called  quadrantal  versors 

Eor  a  long  time  Sir  William  Hamilton  was  not  able  to  find  an  appropriate 
way  to  multiply  quaternions.  While  walking  along  a  canal  in  Dublin  he  discovered  the 
simple  approach,  which  will  be  stated  below,  and  with  his  excitement  he  carved  the 
formula  in  the  stone  of  a  bridge  crossing  the  canal.  Today  a  plaque  marks  this  spot. 

Since 

ij  =k,jk=i  ,ki=j  ; 

Therefore; 

(ij)k  =  kk  =  k2  =  -l, 

(jk)i  =  ii  =  i^  =  -l, 

(ki)j=jj=f=-l 
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as  also 


iOk)  =  ii  =  i^  =  -l, 
j'(ki)=jj=f  =  -l, 
k(ij)  =  kk  =  k^=-l. 

We  therefore  can  omit  the  parenthesis  and  write 


ijk  =  jki  =  kij  =  -l  (Eq.6) 

c.  Representation  of  Quaternions  by  Quadrantal  Versors 

By  using  the  results  in  Equation  6  and  Table  1,  quaternions  can  be  denoted 
as  quadrinomial  expressions,  of  which  one  term  is  called  the  real  part,  while  the  three 
other  terms  made  up  together  a  trinomial,  which  was  called  the  imaginary  part  of  the 
quaternion:  the  square  of  the  former  part  (or  term)  being  always  a  positive,  but  the  square 
of  the  latter  part  (or  trinomial)  being  always  a  negative  quantity.  More  particularly,  this 
imaginary  trinomial  was  of  the  form  ix  +  jy  +  kz,  in  which  x,  y,  z  were  three  real  and 
independent  coefficients,  or  constituents,  and  were,  in  several  applications  of  the  theory, 
constructed  or  represented  by  three  rectangular  coordinates;  while  i,  j,  k  were  certain 
imaginary  units,  or  symbols,  subject  to  the  following  laws  of  combination  as  regards  their 
squares  and  products. [WIEKINS99] 

d.  Identical  Quaternions 

Assume  qi  and  qi  are  two  quaternions. 


qi  =  q2,  (Eq.7) 


in  which 

qi  =  wi  +  ixi  +jyi  +  kzi , 

(Eq.8) 

and 

q2  =  W2  +  ix2  +  jy2  +  kz2 

(Eq.9) 
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Implies 

Wi  =  W2  xi  =  X2  yi  =  y2  and  z\=  Z2 . 

e.  Addition  and  Subtraction  of  Quaternions 

Quaternions  are  added  or  subtraeted  by  adding  or  subtracting  their 
constituents,  so  that 

qi  +  q2  =  (wi  +  W2)  +  i(xi  +  X2)  +  j'(yi  +y2)  +  k(zi  +  Z2)  (Eq.lO) 

f  Multiplication  of  Quaternions 

Multiplication  for  the  primitive  elements  i,  j,  and  k  is  defined  in  Equation 

6. 

Multiplication  of  quaternions  is  defined  by 

qoqi  =  (wO  +  xOi  +  yOj  +  zOk)(wl  +  xli  +  ylj  +  zlk) 

=  (wowi  -  xoXi  -  yoyi  -  zqZi) 

(woxi  +  xowi  +  yozi  -  zoyi)i+ 

(woyi  -  XoZi  +  yoWi  +  zoXi)j+ 

(wozi  +  xoyi  -  yoxi  +  zowQk.  (Eq.  1 1) 

Multiplication  is  not  commutative  in  that  the  products  qoqi  and  qiqo  are 
not  necessarily  equal.  [MORSE53] 

A  simple  but  useful  function  is  the  selection  function 

W(q)  =  W(w  +  xi  +  yj  +  zk)  =  w 

which  selects  the  “real  part”  of  the  quaternion. 

The  quaternion  q  =  w  +xi  +yj  +zk  may  also  be  viewed  as 

q  =  w  +  v 
V  =  xi  +  j'y  +  kz 
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If  we  identify  ^  with  the  3D  vector  (x,  y,  z),  then  quaternion 
multiplication  can  be  written  using  vector  dot  product  ( •  )  and  cross  product  (  x  )  as 

{Wq  +  Vo)(>^i  +  i^i)  =  -  Vq  •  Vj)  +  WqVj  +  WjVq  +  Vq  x  Vj  (Eq.l2) 

In  this  form  it  is  clear  that  qoqi  =  qiqo  if  and  only  if  Vq  x  Vj  =  0 
(i.e.  vectors  are  parallel).  [EBERLY99] 

A  quaternion  q  may  also  be  viewed  as  a  4D  vector  (w,  x,  y,  z).  The  dot 
product  of  two  quaternions  is 

qo  •  qi  =  wowi  +  xqXi  +  yoyi  +  zqZi  =  W(qoqi*  ).  (Eq.l3) 

g.  Conjugate  and  Norm  of  a  Quaternion 

The  conjugate  of  a  quaternion  is  defined  by 

q*  =  (w  +  xi  +  yj  +  zk)*  =  w  -  xi  -  yj  -  zk.  (Eq.  14) 

Therefore  (p*)  *=  p 

The  conjugate  of  a  product  of  quaternions  satisfies  the  properties 

*  *  * 

(pq)  =  q  p  . 

The  norm  of  a  quaternion  is  defined  by 

N(q)  =  N(w  +  xi  +  yj  +  zk)  =  w^  +  x^  +  y^  +  z^.  (Eq.  15) 

Therefore  N(q*)  =  N(q) 

The  norm  is  a  real-valued  function  and  the  norm  of  a  product  of 
quaternions  satisfies  the  properties 

N(pq)  =  N(p)N(q). 

h.  Multiplicative  Inverse  of  a  Quaternion 
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The  multiplicative  inverse  of  a  quaternion  q  is  denoted  q'  ^  and  has  the 
property  qq'  ^  =  q'  ^  q  =  1 .  It  is  construeted  as 

q-'  =  q  /N(q)  (Eq.l6) 

Where  the  division  of  a  quaternion  by  a  real-valued  sealar  is  just 
eomponent  wise  division.  The  inverse  operation  satisfies  the  properties 

(p-  ^  )■ '  =  p  and  (pq)'  ^  =  q  ^  p'  V  (Eq.  17) 

L  Unit  Quaternion  and  Exponential  Form  of  a  Quaternion 

A  unit  quaternion  is  a  quaternion  q  for  whieh  N(q)  =  1 .  The  inverse  of  a 
unit  quaternion  and  the  product  of  unit  quaternions  are  themselves  unit  quaternions.  A 
unit  quaternion  can  be  represented  by 

q=cos^  +  Msin^  (Eq.l8) 

where  u  as  a  3D  vector  has  length  1 .  However,  observe  that  the  quaternion  product 

u  u  -  -  1 .  Note  the  similarity  to  unit  length  complex  numbers 

cos^  +  sin^  .  In  fact,  Euler’s  identity  for  complex  numbers  generalizes  to 

quaternions, 

exp(M^)  =  cos^  +  Msin^  (Eq.l9) 

where  the  exponential  on  the  left-hand  side  is  evaluated  by  symbolically  substituting  ^  ^ 
into  the  power  series  representation  for  exp(x)  and  replacing  products  u  u  by  -  Erom 
this  identity  it  is  possible  to  define  the  power  of  a  unit  quaternion, 

q  *  =  ( cos  ^  +  M  sin  ^  =  QXTp{ut0)  =  cos  t6  +  u  sin  tO  .  (Eq.20) 

It  is  also  possible  to  define  the  logarithm  of  a  unit  quaternion, 

log(q)  =  log( cos ^  +  M  sin ^ )  =  log( exp(M 6)  )=u6  .  (Eq.2 1 ) 

It  is  important  to  note  that  the  noncommutativity  of  quaternion 
multiplication  disallows  the  standard  identities  for  exponential  and  logarithm  functions. 
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The  quaternions  exp(p)  exp(q)  and  exp(p  +  q)  are  not  necessarily  equal.  The  quaternions 
log(pq)  and  log(p)  +  log(q)  are  not  necessarily  equal. 

3,  Application  of  Quaternions  to  Geometric  Rotations 

A  unit  quaternion  q  =  cos  6  +  u  sm  6  represents  the  rotation  of  the  3D  vector 

V  by  an  angle  2  6  about  the  3D  axis  u  .  The  rotated  vector,  represented  as  a  quaternion,  is 
R(  V  )  =  qv  q*  .The  proof  requires  showing  that  R(v  )  is  a  3D  vector,  a  length-preserving 
function  of  3D  vectors,  a  linear  transformation,  and  does  not  have  a  reflection 
component. 

To  see  that  R(  v )  is  a  3D  vector, 

W(R(f  ))=  W(qf  q*) 

=  [(qvq*)  +  (qvq*  )*]/2 
=  [qv  q*  +  qv  *q*  ]/2 
=  q[(v  +  vV2]q* 

=  qW(v)q* 

=  W(T) 

=  0. 

To  see  that  R(  v  )  is  length-preserving, 

N(R(v))  =  N(qvq) 

=  N(q)N(f)N(q*) 

=  N(q)N(f)N(q*) 

=  N(T). 

To  see  that  R(v)  is  a  linear  transformation,  let  a  be  a  real-valued  scalar  and  let 

V  and  ivbe  3D  vectors;  then 

R(a  V  +  >v )  =  q(a  v  +  >v )  q* 
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=  (qavq*)  +  (q  w  q*) 

=  a(qv  q*)  +  (q  w  q*) 

=  aR(v  )  +  R(h’), 

thereby  showing  that  the  transform  of  a  linear  eombination  of  vectors  is  the  linear 
combination  of  the  transforms.  The  previous  three  properties  show  that  R(v)  is  an 
orthonormal  transformation.  Such  transformations  include  rotations  and  reflections. 
Consider  R  as  a  function  of  q  for  a  fixed  vector  v  .  That  is,  R(q)  =  qv  q*. 

This  function  is  a  continuous  function  of  q.  For  each  q  it  is  a  linear  transformation 
with  determinant  D(q),  so  the  determinant  itself  is  a  continuous  function  of  q.  Thus, 

liniq^o  R(q)  =  / ,  the  identity  function 

(the  limit  is  taken  along  any  path  of  quaternions  which  approach  the  zero 
quaternion) 

liniq^o  Z)(q)  =  Z)(0)=  1. 

By  continuity,  D(q)  is  identically  1  and  R(q)  does  not  have  a  reflection 
component.  Now  we  prove  that  the  unit  rotation  axis  is  the  3D  vector  u  and  the  rotation 
angle  is  2  ^  .  To  see  that  m  is  a  unit  rotation  axis  we  need  only  show  that  u  is  unchanged 
by  the  rotation.  Recall  that  u  =  u  u  =  -  1,  which  implies  that  u  =  -  u  . 

R(u)  =  qu  q* 

=  (cos  0  +  iisin  0)u  (cos0-  usin0) 

=  (cos  O  fu  -  (sin  O 
=  (cos  Of'u-  (sin  0f{-  u  ) 

=  u . 

To  see  that  the  rotation  angle  is  2^ ,  let  u  ,  v  ,  and  w  be  a  right-handed  set  of 
orthonormal  vectors.  That  is,  the  vectors  are  all  of  unit  lengthy  •v=y»iv  =  v»w=0; 
and  u  X  V  =w ,  V  X  w  =  u  ,  and  w  x  u=  v  .  The  vector  v  is  rotated  by  an  angle  (j)  to 
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the  vector  qv  q*,  so  v  •  (qv  q*)  =  cos{^).  Using  equation  (13)  and  v  =  -  v  ,  and  =  -  I 
for  unit  quaternions  with  zero  real  part, 

cos(^zi)  =  V  •  (qv  q*) 

=  W(v  *qv  q*) 

=  W[  -  V  (cos  0  +  usm0)v  (cos0  -  iisin^)] 

=  W[(  -  V  cos^  -  V  u  sin0  )(v  cos  0  -v  usin  0  )\ 

2  2  2  2  2 

=  W[  -V  (cos  0  )  +  V  u  sin^  cos  0  -  V  u  V  sin^  cos0  +  (v  u  )  (sin^  )  ] 

=  W[(cos  0  f'-  (sin  ^  )^-  ( ii  +  V  ii  V  )  sin  ^  cos  0  ] 

Now  V  u  =  -  V  •  u+v  xu  =  V  xu  =  -w  and  v  u  v  =  -  w  v  =w  •  v  -w  x  v  =  u  . 
Consequently, 

2  2 

cos(  )  =  W[(cos  ^  )  -(sin^)  -(  u  +  v  u  V  )  sin^  cos  0  ] 

=  W[(cos  0  )^-  (sin  0  )^-  u  (2  sin  0  cos  0  )] 

=  (cos  0  )^-  (sin  0  )^=  cos(2  0  ) 
and  the  rotation  angle  is  (j)  =10  . 

It  is  important  to  note  that  the  quaternions  q  and  -  q  represent  the  same  rotation 
since  ( -  q)  v  ( -  q)*  =  qv  q*.[SHOEMAKE87] 

While  either  quaternion  will  do,  the  interpolation  methods  require  choosing  one 
over  the  other. 

B,  SOURCELESS  MARG  SENSOR  BODY  TRACKING  SYSTEM 

1.  Current  Body  Tracking  Systems 

In  general,  position  and  orientation  tracking  has  seen  inadequate  innovation  and 
development  over  the  past  decade.  This  continues  to  obstruct  advanced  development  of 
immersive  systems  that  allow  participants  to  enter  and  navigate  simulated  environments. 
Today’s  commercial  motion  tracking  systems  are  based  on  optical,  magnetic  and  acoustic 
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sources.  Inertial  sensing  has  been  used  for  head  tracking.  RF  positioning  shows  promise, 
but  no  small-scale  eommercial  systems  are  currently  available  for  indoor  use.  The  most 
popular  traekers  are  aetive  AC  or  DC  magnetic  systems.  [BACHMANN2000]. 

a.  Mechanical  Trackers 

Mechanieal  tracking  systems  are  perhaps  the  oldest  motion  traeking 
technology.  They  offer  the  best  means  of  providing  haptie  feedbaek  to  the  user  of  a 
virtual  environment.  These  systems  are  fairly  aeeurate  and  have  low  latency.  Current 
researeh  generally  involves  using  these  tracking  systems  to  calibrate  other  types  of 
trackers.  [BACHMANN2000] 

Most  famous  example  of  meehanical  body  tracking  system  is  Exoskeleton. 

b.  Optical  Trackers 

Optical  trackers  are  based  on  traeking  of  lumineseent  objeet  by  speeial 
optie  sensors.  (Fig. 3)  This  teehnology  is  known  as  Instant  Marker  Reeognition 
Teehnology  (IMR).  These  systems  usually  consist  of  a  cage,  IR-sensors,  and  data 
processing  computer.  There  are  hundreds  of  digital  sensors  and  eleetronics,  whieh  are 
embedded  in  the  bars  of  cage.  They  employ  IMR  technology  to  cleanly  track  markers. 
Markers  send  optical  signals  for  different  time  intervals.  Since  the  cameras  are  located  in 
fixed  modular  bars,  a  blocked  marker  is  instantly  reacquired  once  occlusion  sends.  This 
produces  consistently  clean  data  streams  with  minimal  data  loss.  (Fig. 3) 


Figure  3.  Optical  Markers  (Courtesy  of  Ascension-teeh  inc  2001) 


17 


Figure  4.  Optical  Motion  Capture  Cage.  Courtesy  of  Ascension-tech  Inc.  2001 

The  disadvantages  of  optical  motion  capture  systems  are  occlusion  of 
body  parts  and  limited  motion  capture  area  (locality). 

c.  Magnetic  Trackers 

Magnetic  trackers  measure  the  position  and  orientation  of  one  or  more 
receiving  antenna  sensors,  typically  located  on  a  user's  head,  hand  or  body,  with  respect 
to  a  transmitting  antenna,  which  is  fixed  in  space.  The  transmitting  antenna  is  driven  by  a 
pulsed,  direct  current  (DC)  signal. 


The  receiving  antenna  measures  not  only  the 
transmitted  magnetic  field  pulse  but  also  the  earth's 
magnetic  field.  A  microprocessor  is  used  to  control  the 
transmitting  and  receiving  elements  and  convert  the 
received  signals  into  position  and  orientation  outputs. 
(Fig.5) 


Figure  5.  Motion  Capture  with  Magnetic  Trackers 
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Electromagnetic  trackers  suffer  from  several  sources  of  error. 
Electromagnetic  interference  (EMI)  from  devices  such  as  radios  or  display  units  can 
cause  erroneous  readings.  Earge  objects  made  of  ferrous  metals  can  interfere  with  the 
electromagnetic  field,  again  causing  inaccuracies  [HAND93]. 

Robustness  is  adversely  effected  by  sensitivity  to  ferromagnetic  objects  in 
the  vicinity,  with  alternating  current  based  trackers  being  more  susceptible  than  direct 
current  based  trackers.  Alternating  current  systems  tend  to  generate  eddy  currents  in 
metallic  objects,  which  then  cause  their  own  electromagnetic  interference  [SKOP96]. 

d.  Fiberoptic  Cable  Deformation  Based  Motion  Capture 

This  system  is  based  on  light  refraction  measurement  techniques.  When 
the  light  travels  through  Eiberoptic  cable  it  has  a  constant  energy  (frequency).  But 
deformation  of  cable  causes  reflection  and  refraction  of  this  beam  inside  the  cable.  And 
one  can  measure  this  energy  changes  and  correlate  with  the  deformation  angle. 

ShapeTape  ©  is  patented  technology  which  uses  this  principle. (Eig. 6) 

ShapeTape  is  a  light  weight,  wearable,  flexible 
ribbon  that  uses  its  own  software,  ShapeWare,  to  create  a  3D 
computer  image  and  data  set  of  its  shape  in  real  time,  based  on 
bend  and  twist  information  from  an  array  of  fiber  optic  sensors 
along  its  length.  It  follows  human  arm,  leg,  back,  and  neck 
movements  for  motion  capture,  virtual  reality,  biomedical, 
gaming,  and  robotic  control  applications.  It  is  also  used  for 
crash  testing,  computer-aided  design,  and  automotive  interior 
design 

Eigure  6.  Shape  Tape© 
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e.  Acoustic  Trackers 

Acoustic  or  ultrasonic  trackers  are  an  inexpensive  option  to  magnetic 
trackers.  They  offer  modest  aecuracies  and  update  rates. 

The  physies  of  sound  limit  the  aeeuracy,  update  rate  and  range  of  aeoustie 
tracking  systems.  Ranges  are  longer  than  that  of  magnetic  trackers  and  magnetie 
interference  is  not  a  problem.  However,  a  clear  line  of  sight  must  be  maintained.  Thus, 
obstruetion  and  shadowing  can  present  difficulties.  Lateney  varies  with  distance  due  to 
the  relatively  slow  speed  of  sound.  Most  eurrent  systems  utilize  40  kHz  tone  pulses. 
Sound  in  this  frequency  band  can  be  rigorously  affected  by  noise  from  metallie  objeets 
sueh  as  key  ehain.  Shorter  wavelengths  more  preeisely  resolve  distanees,  but  quiekly 
attenuate.  In  addition,  high  frequency  omnidirectional  radiators  are  expensive  to 
implement  and  require  more  power. 

f.  Pattern  Recognition 

Passive  stereo  vision  systems  employ  ambient  light  and  square-array 
eharge-eoupled  deviee  (CCD)  eameras  [DURL95].  Multiple  images  from  cameras  with 
varying  viewpoints  are  eompared.  Triangulation  is  then  used  to  determine  position 
[SKOP96].  An  essential  issue  is  to  solve  the  eorrespondenee  problem  of  relating  the  same 
points  in  two  different  images.  Passive  stereo  vision  systems  are  unlikely  to  be  useful  in 
virtual  environment  in  the  near  term,  as  robustness  and  aeeuraey  are  not  yet  eomparable 
to  active  ranging  systems  [DURL95]. 

2.  Fundamental  Principles  of  MARG  Body  Tracking  System 

For  real-time  human  limb  segment  traeking,  most  eurrent  methods  use  aetive 
sources  of  some  sort  to  traek  either  angles  or  the  position  of  referenee  points  on  limbs. 
While  this  has  worked  well  for  some  applieations  sueh  approaehes  typieally  either 
seriously  impede  the  objeets  or  work  only  for  loeal  areas. [YUN99]. 

Reeent  advances  in  miero  maehine  and  embedded  electronie  teehnology  have 
made  it  possible  to  design  a  sourceless  sensor  system.  This  system  uses  earth  magnetie 
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and  gravitational  field  as  references  in  order  to  determine  the  orientation  of  an  object.  An 
example  of  this  system  is  explained  in  detail  for  unmanned  submarine  navigation  in 
[YUN99]. 

MARG  Body  Tracking  system  uses  the  same  idea  with  extended  sensor 
collection.  They  are  magnetic  and  gravitational  field  of  earth  and  angular  rates  of  body 
segments. 

a.  Quaternion  Filter 

MARG  Body  tracking  system  uses  quaternion  angle  sets  to  determine  the 
orientation.  Because  quaternions  are  free  from  singularity  problems  of  Euler  angles  and 
much  more  cheaper  for  rate  transformations. [MCGHEE99] 

Eet  p,  q,  r  be  roll  rate,  pitch  rate  and  yaw  rate  respectively  [COOKE92], 
and  let  q  =  (0  pi  qj  rk)  Then 


e  =  te®fe)  (Eq.22) 

We  can  have  an  exact  orientation  by  integrating  equation  22.  But  such  an 
accurate  and  small  rate  sensor  is  not  available  for  body  tracking  in  today’s  technology. 
Instead,  by  using  micro  machine  technology,  vibrating  beam  sensors  are  more  suitable  for 
this  purpose.  Because,  they  are  small,  cheap,  and  require  negligible  electric  power.  But 
their  error  rates  become  unacceptable  after  only  few  seconds  of  integration.  Eor 
preventing  this  [BACHMANN99]  shows  a  system  with  drift  correction.  Eigure  7  is  the 
quaternion  filter,  which  is  used  in  our  MARG  body  tracking  system. 
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Figure  7.  Quaternion  based  orientation  filter. 

Drift  eorreetion  is  provided  by  Aeeelerometers  and  Magnetometers. 
[YUN2000]  explained  how  to  find  modeling  error  and  estimate  quaternion  orientation  by 
non-linear  regression  analysis  in  detail. 

C.  VIRTUAL  HUMANS  IN  VIRTUAL  ENVIRONMENTS 

1.  Purpose  of  a  Virtual  Human  in  Cyberspace 

Old  religious  seripts  say  that  whole  universe  is  ereated  for  humans.  New 
seriptures  say  that  whole  Cyberspaee  is  ereated  for  Virtual  Humans.  The  main  purpose  of 
a  virtual  human  is  to  represent  human  existenee  within  the  virtual  world. 

The  main  reason  for  developing  realistie  Virtual  Human  is  to  be  able  to  use  them 
in  virtually  any  seene  [THALMANN97].  Potential  human  aetivities,  whieh  may  involve 
humanoids,  are  as  follows 
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•  Simulation  based  learning  and  training 

•  Simulation  of  ergonomie  work  environments 

•  Virtual  patient  for  surgery 

•  Orthopedy  and  prostheses  and  rehabilitation 

•  Architectural  simulation  with  people 

•  Computer  games  involving  people 

•  Game  and  sport  simulation 

•  Interactive  drama  titles  (  play  with  interacting  virtual  characters) 

•  Military  training  simulations 

•  Creating  virtual  terrorist  agents  with  unpredictable  reaction  in  order  to 
prepare  safety  measures. 

•  Crowd  simulations 

•  Virtual  meetings 

Virtual  Humans  will  represent  the  people  in  different  virtual  environments. 
Telepresence  is  the  future  of  multimedia  systems.  This  will  allow  people  to  share 
experiences  and  interact  with  each  other  via  virtual  humans.  These  types  of  virtual 
environments  may  contain  human  controlled  virtual  humans  or  autonomous  humanoids 
or  partially  intelligent  virtual  characters. 

There  are  different  aspects  of  virtual  human  research.  Following  table  shows 
some  of  the  ongoing  projects  within  the  computer  graphics  community. 
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Face  and  body  representation 
Motion  eontrol 


Interaction  with  objects 


Lip  synchronization 


Autonomous  behavior 


Deformable  body  shape 


Interaetion  between  virtual  humans 


Networked  virtual  Environments 


Crowds 


Rendering 


Standard  data  protoeols  for  body  orientations 


Interaetion  by  users 


Optimized  rendering 


Table  2.  Aspeets  of  Virtual  Human  Researeh  [From  THALMANN97] 

Our  research  deals  with  construction  of  virtual  human  with  deformable 
skin  for  a  motion  eapture  system. 

2.  Current  Techniques  For  Inserting  Virtual  Humans  Into  Virtual 
Environments 

Current  techniques  for  inserting  humanoids  into  virtual  worlds  are  not  real-time. 
They  usually  have  pre-produetion  and  postproduetion  phases.  Some  of  these  teehniques 
involve  user  interaetion  to  avatars  by  some  kind  of  device  such  as  mouse  or  joystick.  The 
main  problem  of  these  systems  is  that  they  don’t  represent  a  realistie  human  motion. 
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People  are  skilled  at  perceiving  the  subtle  details  of  human  motion.  A  person  is 
often  able  to  recognize  other  people  at  a  distance  by  their  walk.  Because  of  this  ability, 
people  have  high  standards  for  animations  that  feature  humans.  For  realistic  computer¬ 
generated  motion,  the  virtual  actors  must  move  with  a  natural-looking  style. 

Today  computer  generated  human  motion  is  needed  for  such  applications  as 
animation,  virtual  environments  and  video  games.  For  example,  trainers  could  use  virtual 
sportsmen  to  motivate  and  teach  athletes.  Video  game  designers  could  create  products 
with  highly  interactive,  engaging  characters.  In  the  Navy,  the  movements  of  a  fighting 
soldier  in  simulated  or  real-time  environment  can  be  tracked.  The  ability  to  simulate 
human  motion  also  has  significant  scientific  applications  in  ergonomics,  and  physical 
rehabilitation.  [HODGINGS2000] 

Although  the  applications  for  synthetic  human  motion  are  many,  specifying 
movement  to  a  computer  is  surprisingly  hard.  Even  a  simple  bouncing  ball  can  be 
difficult  to  animate  convincingly,  in  part  because  people  quickly  pick  out  action  that  is 
unnatural  or  implausible  without  necessarily  knowing  exactly  what  is  wrong.  Animation 
of  a  human  is  especially  time-consuming  because  subtleties  of  the  motion  must  be 
captured  to  convey  personality  and  mood. 

There  are  there  types  of  human  motion  generation: 

•  Keyframing 

•  Motion  Capture 

•  Simulation 

Keyframing  forces  the  animator  to  smooth  between  frames.  Simulations  generate 
motion  in  a  fairly  automatic  fashion  but  strictly  depend  on  minimization  of  error  between 
the  physical  based  model  of  motion  and  real  motion  by  computers.  However  body 
tracking  is  a  discrete  sampling  of  real-time  motion  and  gets  as  realistic  as  sampling 
frequency  increases[HODGINS2000] . 
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III.  MODELING  A  VIRTUAL  HUMAN 


A.  KINEMATIC  MODELS 

1.  Introduction 

The  virtual  human  body  is  typically  modeled  using  a  hierarchical  structure.  Joints 
are  the  nodes  of  this  structure  and  segments  such  as  arms,  legs  are  links  between  these 
joints.  In  order  to  create  a  realistic  human  representation  one  needs  to  define  kinematic 
properties  of  this  structure. 

Kinematic  models  specify  geometric  and  time  related  properties  of  the  links 
between  joints.  Position,  orientation,  velocity  and  acceleration  are  examples  of  properties 
of  motion  of  each  link.  Rigid  body  transformation  is  a  fundamental  approach  to  calculate 
orientation  of  each  segment.  End-effector  is  the  final  segment  of  hierarchical  structure. 

There  are  two  types  of  Kinematics.  Forward  Kinematics  determines  the 
orientation  of  end-effector  by  using  the  properties  of  each  segment.  Inverse  kinematics 
determines  the  properties  of  each  segment  in  order  to  reach  final  transformation  of  end- 
effector. 

This  research  used  H-Anim  specification  to  define  kinematic  properties  of  a 
virtual  human.  H-Anim  is  a  structural  model  with  segment  and  joint  containers. 

2,  Rigid  Body  Orientation 

Human  posture  can  be  modeled  as  group  of  orientations  of  body  segments.  Each 
segment  is  a  rigid  body  in  terms  of  physics.  There  are  different  ways  to  represent  rigid 
body  rotations.  This  part  of  this  chapter  will  explain  two  of  them:  Euler  and  quaternion 
rotations. 

a.  Euler  Angle  Rotations 

One  of  the  popular  methods  for  representing  orientation  in  Earth 
coordinates  is  the  Euler  Method.  Euler  method  uses  simple  sequence  of  three  angles.  This 
makes  Euler  Method  an  intuitive  description.  Although  eleven  other  possibilities  exist, 
these  angles  typically  consist  of  the  familiar  azimuth  angle,  y/ ,  the  elevation  angle,  6  , 
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and  the  roll  angle,  (p  ,  [CRAI89].  Euler  angles  speeify  suecessive  rotations  to  bring  the 
Earth  eoordinates  into  alignment  with  the  body  eoordinates  [COOK92]. 

Euler’s  theorem  states  that  any  numbers  of  rotations  of  a  rigid  body  about 
an  Earth-fixed  axis  are  equivalent  to  a  single  rotation  about  a  single  Earth-  fixed  axis.  If 
all  rotations  are  about  the  north(x),  east(y)  and  down(z)  axes,  as  depleted  in  Eigure  8  ,  the 
angles  are  e ailed  Euler  Angles. 


Eigure  8.  Right  handed  Euler  Rotations  (Earth  Axis) 


The  standard  “azimuth,  elevation,  and  roll”  set  is  represented  by  the  three 
rotation  matrices  given  in  Equation  23.  Rx  represents  a  rotation  about  the  x  axis 
[CRAI89],  and  similarly  for  Ry  and  Rz. 

R  =  RzRyRx  (Eq.23) 

These  symbols  are  reserved  for  Euler  Rotations  with  the  angles  \i/,0,(j) . 
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Right  Hand  rule  determines  the  direction  (sign)  of  the  rotation.  Limits  of  the  rotations  are 
as  follows: 


y/  =  ±7r. 


.=4 


(j)  =  ±7r. 


(Eq.24) 


X-axis  rotation  (roll)  matrix  is  given  by 

"10  o' 

R^{(p)=  0  cos^3  sin^  (Eq.25) 

0  -sin^  cos^3 

Y-axis  rotation  (elevation,  pitch)  ; 

cos^  0  -sin  6* 

R^{e)=  0  1  0  (Eq.26) 

sin  6*  0  cos  6* 

Z-axis  rotation  (azimuth,  heading,  yaw)  is  given  by 

cos^i^  sin^  0 

R^{y/)=  -sin^  cos^  0  (Eq.27) 

0  0  1 

Every  rotation  matrix  is  a  3x3  matrix,  but  to  make  both  translation  and 
rotation  calculations  using  matrices,  4x4  homogeneous  matrices  are  used.  A 
homogeneous  rotation  matrix  in  graphics  takes  the  form  [EOLE97] 


h  r,  a 

r,  r,  rg  b 

r,  rg  c 

0  0  0  1 


(Eq.28) 
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Rotation  of  a  vertex  is  calculated  by  multiplying  these  rotation  matrices 
with  homogeneous  coordinates  of  this  vertex.  First  rotation  is  associated  with  the 
rightmost  matrix.  And  successive  multiplications  determine  the  final  orientation  of  the 
object.  For  example  a  roll  followed  by  an  elevation  and  another  azimuth  rotation  will 
form  matrix  R  by  multiplying  corresponding  matrices  from  right  to  left. 

R  =  RzRyRx  (Eq.29) 


R  = 


cos  cos  6* 
sin  ^  cos  6* 
-  sin  6* 


sin  (p  sin  6  cos  y/  -  cos  cp  sin  y/ 
sin  (p  sin  6*  sin  ^  +  cos  (p  cos  y/ 
cos  6*  sin  ^ 


sin  ^  sin  ^  +  cos  cp  cos  y/  sin  6 
cos  (p  sin  ^  sin  ^  -  cos  y/  sin  cp 
cos  ^3  cos  ^ 


(Eq.30) 


b.  Body  Rates  to  Euler  Rates 

Representing  a  rotational  velocity  around  the  nose  vector,  side  vector(s), 
and  the  approach  or  belly  vector  requires  a  standard  frame.(Eig.9)  Two  common 
coordinate  frames  are  Earth  and  Body  coordinate  frames. 


Rotational  velocity  in  Earth  coordinates  is; 


(Eq.31) 
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Figure  9.  Body  Coordinates  Axis 

And  in  body  coordinates 

Bco=[p  q  rj  (Eq.32) 

where  p  is  roll  rate,  q  is  pitch  rate  and  r  is  yaw  rate.  These  symbols  are  reserved  for  body 
rotations  [MCGHEE93].  Euler  rates  are  in  Earth  coordinates  while  p.q.r  are  in  body 
coordinates. 

After  doing  calculation  to  solve  Body  coordinates  in  terms  of  Earth 
coordinates  [DUMAN2000]  we  have  the  following  results. 

(p  1  tan  6*  sin  ^  tan  6*  cos  ^  p 

6  =  Q  cos(p  -sin^  q  (Eq.33) 

\(/  0  sec  6*  sin  ^  sec  6*  cos  ^3  r 

Equation  33  is  known  as  the  gimbal  rate  equations  and  the  matrix  is 
named  as  T  matrix.  Currently  this  is  the  most  used  form  for  animation  and 
simulations.  [ERAN69] 
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c.  Rigid  Body  Dynamics  with  Quaternions 

Euler  angles  represent  the  rotation  of  rigid  body  around  three  orthogonal 
coordinates.  However  we  can  describe  the  orientation  of  these  objects  by  a  single  rotation 
{6)  around  a  single  vector( u ).  Moreover  we  can  constraint  this  vector  with  a  unit 
magnitude.  Thus  unit  quaternion  representing  this  rotation  is  [MCGHEE97] 


follows, 


9  = 


e 


.  e 


cos—  MSin  — 


(Eq.34) 


2  Ij 

[SHOEMAKERS?]  defines  the  rotation  of  vector  p  by  quaternion  q  as 


rotated 


=  qpq 


(Eq.35) 


The  reason  of  half  angle  in  Equation  34  is  that  quaternion  q  in  Equation  35 
rotates  the  vectors  perpendicular  component  twice  the  angle  which  is  perpendicular  to 
rotation  axis. 

After  redefining  quaternion  rotations  now  we  need  to  derive  quaternions 
from  body  rates.  Angular  rates  p  ,q,  r  can  be  used  to  form  a  q  quaternion  rate  relative  to 
Earth  coordinates. 


Assumptions: 


0  ,  .00 

cos— =  1  ,  sm— =  — 

2  2  2 


Thus 


(  ^ 

^  .  0^ 

cos  — 

u  sm  — 

=  q  =  \ 

V  — 

1  2 

2? 

2j 

(Eq.36) 


(Eq.37) 
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For  a  small  time  change,  the  rotation  quaternion  can  be  expressed  as  a 
function  of  time  as  follows. 


f 


q{t)  =  1^1 


1  .  ^ 

-vet 

2  ) 


(Eq.38) 


is  angular  rate  about  a  vector  v  in  body  vectors  which  is  equivalent  to 

body  rates  p,  q,  r. 


and 


F^=(p  qr) 


(Eq.39) 


r  1  1  1  ^ 

q{t)=  1  -pt-qt-rt 
\  2  2  Z  J 


(  111^ 
q{t)=  0  -p-q-r 
\  2  2  2  ) 


4(0  =  ^(0  P  q  r) 


(Eq.40) 


m  =  \B.  (Eq.41) 

Eet  qi  be  initial  orientation  in  Earth  coordinates  and  qi  be  a  second 
rotation  in  body  coordinates  and  q3  be  the  composite  of  these  two 
rotations.  [B  ACHMANN2000] 


^3  =q\q2 
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Let’s  derive  qs  with  respeet  to  time. 


^3  +^2 

Derivation  of  qi  will  be  zero  beeause  it  is  eonstant.  And  by  using  Equation 
41  derivation  of  q2  can  be  substituted  by  its  rate  values. 

^3  (0  =  ^1 1 (Eq.42) 
d.  Rotation  Matrices  From  Quaternions 

Currently  there  isn’t  a  computer  graphics  hardware,  which  is  directly 
using  quaternions  for  displaying  vertices  on  the  screen.  This  forced  virtual  human 
interface  to  own  a  converter  from  quaternion  rotation  to  axis  angle  rotations.  Conversion 
can  be  achieved  by  using  following  matrix. 

Eet  q  =  (w  +  xi  +  yj  +  zk)  then  , 

w^+x^-y^-z^  2xy  +  2wz  2xz-2wy  0 

2xy-2wz  w^-x^+_y^-z^  2xy  +  2wx  0 

2xz  +  2wy  2_yz-2>vx  w^-x^-y^+z^  0 

0  0  0  +  z^ 

Rotated  vector  can  be  calculated  by  multiplication  of  Q  matrix  to  vector 

P. 

3.  Comparison  Of  Quaternion  And  Euler  Angles  For  Kinematic  Models 

Quaternions  and  Euler  angles  each  have  their  own  advantages  and  disadvantages. 
The  most  major  advantage  of  quaternions  is  that  no  singularity  exists  when  the  elevation 
angle  {6)  passes  through  ±;r  12.  In  the  Euler  Method,  roll  and  azimuth  Euler  angle  rates 
become  very  large  (infinite)  due  to  division  by  zero.  Truncating  the  angles  at  ±;z'/2  will 
avoid  this  problem.  However,  this  truncating  will  result  in  jump  during  the  rotation 
[COOK92]. 
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Quaternion  rotation  avoids  the  “braneh  cut”  problem  of  Euler  angles.  When  using 
quaternions  to  rotate  an  object,  after  every  full  rotation  the  quaternion  representing  the 
rotation  will  switch  to  its  negative. 

Quaternions  can  be  computed  directly  from  the  dynamic  equations,  bypassing  the 
computation  of  transcendental  functions  necessary  in  computing  Euler  angles.  This  direct 
arithmetic  operation  reduces  the  cost  of  computation  compared  to  matrix  multiplication. 

Quaternions  are  compact  and  simple.  However,  there  are  three  difficulties  with 
using  quaternions. 

•  Each  orientation  of  an  object  can  actually  be  represented  by  two  quaternions. 
Since  rotation  about  the  axis,  v  by  an  angle  0  is  the  same  as  rotation  about  -v  by 
the  angle  -0  ;  the  corresponding  quaternions  are  antipodal  points  on  the  sphere  in 
4D. 

•  Orientations  and  rotations  are  not  exactly  the  same  thing.  In  an  animation,  360“ 
rotation  is  very  different  from  a  rotation  of  0“.  With  the  first  one,  the  model  will 
be  animated  for  a  full  rotation.  In  the  second  one,  no  animation  will  be  executed. 
After  the  rotation  however,  the  same  quaternion  (1  0  0  0)  or  (-1  0  0  0)  represents 
both.  Thus  specifying  multiple  rotations  with  quaternions  requires  intermediate 
control  points  [EOEE97]. 

•  Quaternions  provide  an  isotropic  method  for  rotation.  The  interpolation  is 
independent  of  everything  except  the  relation  between  the  initial  and  final 
rotations.  This  is  useful  for  interpolating  orientations  of  tumbling  bodies,  but  not 
for  interpolating  the  orientations  of  a  virtual  camera  in  a  scene.  Humans  strongly 
prefer  cameras  to  be  upright,  and  are  profoundly  disturbed  by  tilted  cameras 
[EOEE97].  By  its  very  nature,  the  notion  of  a  preferred  direction  cannot  easily  be 
built  into  the  quaternion  representation  [WATT98].  Quaternions  have  no  such 
preferences,  and  therefore  should  not  be  used  for  camera  interpolation.  The  lack 
of  an  adequate  method  for  interpolating  complex  camera  motion  has  led  to  many 
computer  animations  having  static  cameras  or  very  limited  camera  motion 
[EOEE97]. 
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B,  H-ANIM  2,0  STANDARD  FOR  HUMAN  MODELS 

1.  Introduction 

As  the  3D  Internet  eontinues  to  grow,  there  will  be  an  inereasing  need  to 
charaeterize  human  beings  in  online  virtual  environments.  Aehieving  that  goal  will 
require  the  establishment  of  libraries  of  interehangeable  humanoids,  as  well  as  authoring 
tools  that  make  it  easy  to  ereate  new  humanoids  and  animate  them  in  different  ways. 

H-Anim  group  speeified  design  goals  as  follows: 

Compatibility:  Humanoids  should  work  in  any  VRML97  eompliant  browser. 

Flexibility:  No  assumptions  are  made  about  the  types  of  applieations  that  will  use 
humanoids. 

Simplicity:  The  speeifieation  ean  always  be  extended  later. 

Mueh  of  the  design  is  driven  by  these  three  goals.  The  eompatibility  requirement 
avoids  the  use  of  seripting,  sinee  the  VRML97  speeifieation  does  not  require  any 
partieular  seripting  language  to  be  implemented.  The  flexibility  requirement  makes  the 
spee  fairly  "low-level",  in  that  it  allows  direct  access  to  the  joints  of  the  humanoid's  body 
and  the  vertices  of  the  individual  body  segments.  The  simplicity  requirement  focuses 
specifically  on  humanoids,  instead  of  trying  to  deal  with  arbitrary  articulated  figures. 

The  human  body  consists  of  a  number  of  segments  (such  as  the  forearm,  hand  and 
foot),  which  are  connected  to  each  other  by  joints  (such  as  the  elbow,  wrist  and  ankle).  In 
order  for  an  application  to  animate  a  humanoid,  it  needs  to  obtain  access  to  the  joints  and 
alter  the  joint  angles.  The  application  may  also  need  to  retrieve  information  about  such 
things  as  joint  limits  and  segment  masses. 
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Figure  10.  HANIM  Joint  Hierarchy 
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A  mesh  of  polygons  will  typically  define  each  segment  of  the  body,  and  an 
application  may  need  to  alter  the  locations  of  the  vertices  in  that  mesh.  The  application 
may  also  need  to  obtain  information  about  which  vertices  should  be  treated  as  a  group  for 
the  purpose  of  deformation. 

An  H-Anim  file  contains  a  set  of  Joint  nodes  that  are  arranged  to  form  a 
hierarchy.  Each  Joint  node  can  contain  other  Joint  nodes,  and  may  also  contain  a  Segment 
node,  which  describes  the  body  part  associated  with  that  joint.  Each  Segment  can  also 
have  a  number  of  Site  nodes,  which  define  locations  relative  to  the  segment.  Sites  can  be 
used  for  attaching  clothing  and  jewelry,  and  can  be  used  as  end-effectors  for  inverse 
kinematics  applications.  They  can  also  be  used  to  define  eye-points  and  viewpoint 
locations. 

Each  Segment  node  can  have  a  number  of  Displacer  nodes  that  specify  which 
vertices  within  the  segment  correspond  to  a  particular  feature  or  configuration  of  vertices. 
[HANIM2001] 

2.  Fundamental  Prototypes  For  H-Anim  2,0  Standard 

There  are  five  fundamental  nodes  in  H-anim  2001  spec.  Here  I  explained  only  3 
of  them.  Because  my  implementation  used  only  Humanoid  ,  joint ,  and  segment  nodes. 

a.  Humanoid  Node 

The  Humanoid  node  serves  as  the  overall  container  for  the  Joint,  Segment, 
Site  and  Viewpoint  nodes,  which  define  the  skeleton,  geometry  and  landmarks  of  the 
humanoid  figure.  Additionally,  the  node  provides  a  means  for  defining  information  about 
the  author,  copyright  and  usage  restrictions  of  the  model.  Eater  in  this  document  more 
specific  details  about  the  interface  for  this  node  will  be  described,  but  for  now  we  will 
give  a  general  overview  of  the  entities  involved  and  the  different  ways  in  which  the 
humanoid  figure  can  be  defined. 

The  geometry  of  a  H-Anim  humanoid  figure  can  be  described  in  two 
ways.  The  first  way  is  within  the  scene  graph  of  the  joint  hierarchy,  which  is  described  in 
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the  skeleton  field  of  the  Humanoid  node.  Geometry  (in  the  form  of  Shape  nodes)  defined 
within  Segment  nodes  of  this  joint  hierarehy  describes  the  body  as  separate  geometric 
pieces.  This  method,  while  computationally  efficient,  has  certain  visual  anomalies  (such 
as  seams  or  creases),  which  detract  from  the  appearance  of  the  humanoid  figure.  The 
second  way  of  describing  the  geometry  of  an  H-Anim  humanoid  figure  is  as  a  continuous 
piece  of  geometry,  within  the  skin  field  of  the  Humanoid  node.  For  this  method,  point 
and  normal  vector  data  sets  are  first  defined  in  the  skinCoord  and  skinNormal  fields  (in 
the  form  of  Coordinate  and  Normal  nodes,  respectively).  This  data  is  defined  in  this 
manner  as  to  allow  it  to  be  referenced  by  two  different  entities  within  the  definition  of  the 
Humanoid  node.  The  first  entity  is  an  IndexedFaceSet  node(s),  which  defines  the 
surface  of  the  geometry  of  the  humanoid  figure  within  the  skin  field.  In  most  cases  this 
surface  will  be  a  single  IndexedFaceSet  node,  however  the  surface  can  also  be  defined  as 
multiple  IndexedFaceSet  nodes.  It  is  possible  that  depending  on  the  implementation  of 
the  IndexedFaceSet  nodes  and  the  configuration  of  the  humanoid  figure,  multiple 
IndexedFaceSet  nodes  may  provide  better  performance  by  isolating  the  continuous  mesh 
changes  to  localized  surfaces. 

A  skeletal  hierarchy  is  defined  for  each  H-Anim  humanoid  figure  within 
the  skeleton  field  of  the  Humanoid  node.  This  hierarchical  definition  of  joints  is  present 
even  when  the  geometry  of  the  humanoid  figure  is  not  defined  within  the  skeleton  field, 
as  the  Joint  definitions  are  still  needed  to  describe  the  coordinate  frame  about  which  the 
vertices  of  the  continuous  mesh  are  deformed.  The  H-Anim  specification's  only 
requirement  is  that  this  skeletal  hierarchy  contains  the  HumanoidRoot  Joint.  All  of  the 
other  Joint  nodes  are  optional  and  are  not  required  for  a  humanoid  figure  to  be  H-Anim 
compliant.  Naturally,  most  H-Anim  figures  will  have  many  more  joints  defined  than  this 
required  one. 
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PROTO  Humanoid  [ 


field  SFVecSf  bboxCenter 

00  0 

field  SFVec3f  bboxSize 

-1  -1  -1 

exposedField 

SFVec3f 

center 

0  00 

exposedField 

MFString 

info 

[] 

exposedField 

MFNode 

joints 

[] 

exposedField 

SFString 

name 

ft  ft 

exposedField 

SFRotation  rotation 

0010 

exposedField 

SFVec3f 

scale 

111 

exposedField 

SFRotation  sealeOrientation  OOF 

exposedField 

MFNode 

segments 

[] 

exposedField 

MFNode 

sites 

[] 

exposedField 

MFNode 

skeleton 

[] 

exposedField 

MFNode 

skin 

[] 

exposedField 

SFNode 

skinCoord 

NULL 

exposedField 

SFNode 

skinNormal 

NULL 

exposedField 

SFVec3f 

translation 

000 

exposedField 

SFString 

version 

"2.0" 

exposedField 

MFNode 

viewpoints 

[] 

Figure  1 1 .  Humanoid  Node  Prototype 


The  bboxCenter  and  bboxSize  fields  speeify  a  bounding  box  that  eneloses 
the  ehildren  of  the  Humanoid  node.  This  is  a  hint  that  may  be  used  for  optimization 
purposes.  A  default  bboxSize  value,  (-1,  -1,  -1),  implies  that  the  bounding  box  is  not 
speeified  and,  if  needed,  shall  be  ealeulated  by  the  browser. 

The  eenter  field  speeifies  a  translation  offset  from  the  origin  of  the  loeal 
eoordinate  system  (0,0,0).  The  rotation  field  speeifies  a  rotation  of  the  eoordinate  system 
of  the  humanoid  figure.  The  seale  field  speeifies  a  non-uniform  seale  of  the  humanoid 
figure  eoordinate  system  and  the  seale  values  must  be  greater  than  zero.  The 
sealeOrientation  speeifies  a  rotation  of  the  eoordinate  system  of  the  humanoid  figure 
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before  the  seale  (to  speeify  seales  in  arbitrary  orientations).  The  sealeOrientation  applies 
only  to  the  seale  operation.  The  translation  field  speeifies  a  translation  to  the  eoordinate 
system  of  the  entire  humanoid  figure. 


The  name  field,  whieh  must  be  present,  stores  the  name  of  the  humanoid 
defined  by  the  Humanoid  node. 

The  skeleton  field  eontains  the  HumanoidRoot  Joint  node.  The  Humanoid 
node  is  eonsidered  the  parent  node  of  the  HumanoidRoot  Joint  node  and  defines  a 
eoordinate  spaee  for  the  HumanoidRoot  Joint  node.  This  means  that  the  transformation 
of  the  Humanoid  node  aeeumulates  down  the  seene  graph  defined  within  the  skeleton 
field.  As  was  mentioned  in  the  introduetion  of  this  seetion  hierarehy  of  Joint  nodes  are 
defined  for  eaeh  H-Anim  humanoid  figure  within  the  skeleton  field  of  the  Humanoid 
node  and  hierarehieal  definition  of  joints  is  present  even  when  the  geometry  of  the 
humanoid  figure  is  not  defined  within  the  skeleton  field. 

The  joints  field  eontains  a  list  of  referenees,  one  for  each  Joint  node 
defined  within  the  skeleton  field  hierarchy  of  the  Humanoid  node.  This  field  was  added 
to  the  H-Anim  specification  because  of  a  functional  limitation  with  VRML97  when 
querying  the  system  for  a  list  of  all  the  Joint  nodes  defined.  It  is  possible  that  in  future 
versions  of  the  H-Anim  specification  this  list  will  not  be  necessary,  however  for  now  all 
implementations  of  H-Anim  humanoid  figures  are  required  to  implement  this  feature.  The 
order  in  which  the  joints  are  listed  is  irrelevant,  since  the  names  of  the  joints  are  stored  in 
the  Joint  nodes  themselves. 

The  segments  field  contains  a  list  of  references,  one  for  each  Segment 
node  defined  within  the  skeleton  field  hierarchy  of  the  Humanoid  node.  This  field  was 
added  to  the  H-Anim  specification  because  of  a  functional  limitation  with  VRML97 
when  querying  the  system  for  a  list  of  all  the  Segment  nodes  defined.  It  is  possible  that 
in  future  versions  of  the  H-Anim  specification  this  list  will  not  be  necessary,  however  for 
now  all  implementations  of  H-Anim  humanoid  figures  are  required  to  implement  this 
feature.  The  order  in  which  the  segments  are  listed  is  irrelevant,  since  the  names  of  the 
segments  are  stored  in  the  Segment  nodes  themselves. 
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The  sites  field  contains  a  list  of  references,  one  for  each  Site  node  defined 
within  the  skeleton  field  hierarchy  of  the  Humanoid  node.  This  field  was  added  to  the  H- 
Anim  specification  because  of  a  functional  limitation  with  VRML97  when  querying  the 
system  for  a  list  of  all  the  Site  nodes  defined.  It  is  possible  that  in  future  versions  of  the 
H-Anim  specification  this  list  will  not  be  necessary,  however  for  now  all 
implementations  of  H-Anim  humanoid  figures  are  required  to  implement  this  feature.  The 
order  in  which  the  sites  are  listed  is  irrelevant,  since  the  names  of  the  sites  are  stored  in 
the  Site  nodes  themselves. 

The  viewpoints  field  has  a  different  functionality  and  behavior  than  the 
joints,  segments  and  sites  fields.  The  viewpoints  field  can  contain  zero  or  more 
Viewpoint  node  definitions.  The  Viewpoint  nodes  defined  within  this  field  are  affected 
by  the  transformations  applied  to  the  Humanoid  node,  but  are  not  affected  by  any  of  the 
transformations  performed  on  the  skeletal  hierarchy  defined  within  the  skeleton  field. 
The  order  in  which  the  Viewpoint  nodes  are  listed  will  determine  what  order  they  appear 
in  the  Viewpoint  stack. 

The  skin  field  contains  a  single  or  a  collection  of  IndexedFaceSet  nodes. 
Those  IndexedFaceSet  nodes  reference  the  Coordinate  and  Normal  nodes  defined  within 
the  skinCoord  and  skinNormal  fields,  respectively,  of  the  Humanoid  node.  The 
Coordinate  and  Normal  nodes  are  defined  in  this  manner  so  that  the  data  within  them  can 
be  accessed  by  two  entities.  The  first  entity  is  the  Humanoid  node  itself,  which  has 
internal  mechanisms  that  can  manipulate  the  Coordinate  and  Normal  data  sets  to  create 
the  appropriate  surface  deformations  as  the  Joint  nodes  of  the  humanoid  figure  are 
animated.  The  second  entity  that  uses  these  data  sets  is  the  IndexedFaceSet  nodes 
defined  within  the  skin  field,  which  describes  the  surface  geometry  of  the  humanoid 
figure. 

The  skinCoord  field  contains  a  single  Coordinate  node  definition,  which  is 
used  by  the  internal  mechanisms  of  the  Humanoid  node  to  create  the  appropriate  surface 
deformations  as  well  as  the  IndexedFaceSet  node  which  does  the  actual  rendering  of  the 
surface  geometry. 
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The  skinNormal  field  contains  a  single  Normal  node  definition  which  is 
used  by  the  internal  mechanisms  of  the  Humanoid  node  to  create  the  appropriate  surface 
deformations  as  well  as  the  IndexedFaceSet  node,  which  contains  the  actual  surface 
geometry  that  is  rendered. 

The  Humanoid  node  should  be  given  a  DBF  name  of  "Humanoid",  in 
order  to  make  it  easy  to  access  using  the  External  Authoring  Interface. 

b.  Segment  Node 

The  Segment  node  is  used  describe  the  attributes  of  the  physical  links 
between  the  joints  of  the  humanoid  figure.  Each  body  part  (pelvis,  thigh,  calf,  etc)  of  the 
humanoid  figure  is  represented  by  a  Segment  node.  These  nodes  are  organized  in  their 
skeletal  hierarchy  of  the  humanoid  and  provide  a  container  for  information  which  is 
specific  to  each  segment  of  the  body. 

The  Segment  node  is  a  specialized  grouping  node  that  provides  a  container 
for  nodes  in  it's  children  field.  See  the  VRME97  specification,  specifically  the  4.6.5, 
Grouping  and  children  nodes  section,  for  more  details  on  grouping  nodes.  The  reason  a 
Segment  node  is  considered  a  specialized  grouping  node  is  that  it  can  only  be  defined  as  a 
child  of  an  H-Anim  Joint  node.  To  do  otherwise  would  be  a  violation  of  the  structure  of 
the  H-Anim  Specification. 

It  is  suggested,  but  not  required,  that  all  of  the  body  segments  should  be 
built  in  place.  That  is,  they  should  require  no  translation,  rotation,  or  scaling  to  be 
connected  with  their  neighbors.  In  such  a  scenario  all  the  body's  coordinates  will  share  a 
common  origin,  which  is  that  of  the  humanoid  itself  If  this  proves  difficult  for  an 
authoring  tool  to  implement,  it  is  acceptable  to  use  a  Transform  node  inside  each 
Segment,  or  even  several  Transform  nodes,  in  order  to  position  the  geometry  for  that 
segment  correctly.  The  reason  it  is  best  to  avoid  such  transformations  is  that  it  will 
impact  runtime  performance.  It  is  also  recommended  that  EOD  nodes  are  used  to 
describe  the  geometry  of  Segment  nodes  of  high  complexity. 
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This  section  describes  the  interface  of  the  Segment  node  and  the  fields 
which  allow  the  functionality  of  the  node  to  be  defined  and  accessed. 


PROTO  Segment  [ 

field  SFVec3f  bboxCenter  0  0  0 
field  SFVec3f  bboxSize  -1  -1  -1 


exposedField 

SFVec3f 

centerOfMass  0  0  0 

exposedField 

MFNode 

children  [  ] 

exposedField 

SFNode 

coord  NULL 

exposedField 

MFNode 

displacers  [  ] 

exposedField 

SFFloat 

mass  0 

exposedField 

MFFloat 

moments Ofinertia  [000000000] 

exposedField 

SFString 

name  "" 

eventin  MFNode  addChildren 

eventin  MFNode  removeChildren 

] 


Figure  12.  Segment  Node  Prototype 


The  bboxCenter  and  bboxSize  fields  specify  a  bounding  box  that  encloses 
the  children  of  the  Segment  node.  This  is  a  hint  that  may  be  used  for  optimization 
purposes.  A  default  bboxSize  value,  (-1,  -1,  -1),  implies  that  the  bounding  box  is  not 
specified  and,  if  needed,  shall  be  calculated  by  the  browser. 

The  name  field  must  be  present,  so  that  the  Segment  can  be  identified  at 
runtime.  The  name  field  is  used  for  identifying  the  Segment  nodes  at  runtime.  Each 
Segment  should  have  a  DEE  name  that  matches  the  name  field  for  that  Segment,  but  with 
a  distinguishing  prefix  in  front  of  it.  That  prefix  can  be  anything,  but  must  be  the  same 
for  all  the  Segment  nodes  in  a  particular  humanoid.  The  distinguishing  prefix  is  useful  in 
the  case  of  static  routing  to  the  Segment  nodes  of  multiple  humanoids  in  the  same  file.  If 
only  a  single  humanoid  is  stored  in  a  file,  the  prefix  should  be  "hanim  "  (for  Humanoid 
Animation).  Eor  example,  the  left  thigh  would  have  a  DEE  name  of  "hanim  l  thigh". 
The  DEE  name  is  used  for  static  routing  and  the  name  field,  which  must  be  present,  is 
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used  by  applications  to  identify  the  Segment  node  at  runtime.  All  the  other  fields  are 
optional,  but  must  be  defined  in  the  PROTOtype. 

The  mass  field  is  the  total  mass  of  the  segment. 

The  centerOfMass  field  is  the  location  within  the  segment  of  its  center  of 
mass.  Note  that  a  zero  value  was  chosen  for  the  mass  in  order  to  give  a  clear  indication 
that  no  mass  value  is  available. 

The  momentsOflnertia  field  contains  the  moment  of  inertia  matrix.  The 
first  three  elements  are  the  first  row  of  the  3x3  matrix,  the  next  three  elements  are  the 
second  row,  and  the  final  three  elements  are  the  third  row. 

The  coord  field  is  used  for  Segment  nodes  that  have  deformable  meshes 
and  should  contain  a  Coordinate  node  that  is  USEd  in  the  IndexedFaceSet  for  the 
Segment.  The  Coordinate  node  should  be  given  the  same  name  DEE  name  as  the 
Segment,  but  with  a  "  coords"  appended  (e.g.  "skull  coords"). 

The  displacers  field  stores  the  Displacer  nodes  for  a  particular  Segment 

node. 

c.  Joint  Node 

The  Joint  node  is  used  as  a  building  block  to  describe  the  articulations  of 
the  humanoid  figure.  A  Joint  node,  each  of  which  is  organized  into  a  hierarchy  that 
describes  the  overall  skeleton  of  the  humanoid,  represents  each  articulation  of  the 
humanoid  figure.  This  hierarchy  describes  the  inherent  parent-child  relationship  of  the 
skeleton  as  well  as  providing  a  container  for  information,  which  is  specific  to  each  joint 
of  the  skeleton.  Within  this  description,  we  will  give  more  specific  details  about  the 
interface  for  this  node,  but  for  now  we  will  give  a  general  overview  of  the  entities 
involved  and  the  different  ways  in  which  the  Joint  nodes  interact  with  the  geometry  of  the 
overall  humanoid  figure.  Before  attempting  to  understand  the  operation  of  the  Joint 
node,  the  reader  should  become  familiar  with  the  Humanoid  node  and  its  operation.. 

The  Joint  node  is  a  specialized  grouping  node  that  defines  a  coordinate 
system  for  nodes  in  its  children  field  that  is  relative  to  the  coordinate  systems  of  its  parent 
node.  The  reason  a  Joint  node  is  considered  a  specialized  grouping  node  is  that  it  can 
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only  be  defined  as  a  child  of  another  Joint  node  or  in  the  case  of  the  HumanoidRoot  Joint, 
as  the  first  node  of  the  skeletal  hierarchy  defined  in  the  skeleton  field  of  the  Humanoid 
node.  Additionally,  a  Joint  node  has  two  fields,  which  allow  it  to  manipulate  individual 
vertices  from  the  skinCoord  field  of  the  Humanoid  node.  Incoming  rotation  field  events 
of  the  Joint  node  affect  the  vertices  indicated  by  the  skinCoordlndex  field  by  a  factor, 
which  is  described  by  the  corresponding  values  within  the  skinCoordWeight  field  of  the 
Joint  node.  The  skinCoordWeight  field  contains  a  list  of  floating  point  values,  which 
describe  an  amount  of  "weighting"  that  should  be  used  to  affect  the  appropriate  indices 
(as  indicated  by  the  skinCoordlndex  field)  of  the  skinCoord  field  of  the  Humanoid  node. 
It  is  worth  pointing  out  that  the  skinCoordWeight  and  skinCoordlndex  fields  are  only 
used  with  continuous  mesh  models. 

The  Joint  node  is  also  used  to  store  other  joint-specific  information.  In 
particular,  a  joint  name  is  provided  so  that  applications  can  identify  each  Joint  node  at 
runtime.  The  Joint  node  may  contain  hints  for  inverse-kinematics  systems  that  wish  to 
control  the  H-Anim  figure.  These  hints  include  the  upper  and  lower  joint  limits,  the 
orientation  of  the  joint  limits,  and  a  stiffness/resistance  value.  Note  that  these  limits  are 
not  enforced  by  any  mechanism  within  the  scene  graph  of  the  humanoid,  and  are 
provided  for  information  purposes  only.  Use  of  this  information  and  enforcement  of  the 
joint  limits  is  up  to  the  application. 

Figure  10  describes  the  interface  of  the  Joint  node  and  the  fields  which 
allow  the  functionality  of  the  node  to  be  defined  and  accessed. 
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PROTO  Joint  [ 

exposedField 
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exposedField 

MFNode 
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[] 

exposedField 
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[] 
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exposedField 

SFRotation 
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exposedField 
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111 

exposedField 
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exposedField 

MFInt32 
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exposedField 

MFFloat 
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exposedField 
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exposedField 
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0  00 
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[] 

] 

Figure  13.  Joint  Node  Prototype 


The  center  field  specifies  a  translation  offset  from  the  root  of  the  overall 
humanoid  body  description  and  is  not  intended  to  receive  events  in  most  cases.  Since  all 
of  the  Joint  nodes  have  their  center  values  defined  in  the  same  coordinate  frame,  the 
length  of  any  segment  can  be  determined  by  calculating  the  distance  between  the  parent's 
Joint  center  and  the  child's  Joint  center.  The  exception  will  be  segments  at  the  ends  of  the 
fingers  and  toes,  for  which  the  end-effector  Site  nodes  within  the  Segment  node  must  be 
used. 

The  rotation  field  specifies  a  rotation  of  the  coordinate  system  of  the  Joint 
node,  fields  necessary 

The  skinCoordWeight  field  contains  a  list  of  floating  point  values,  which 
describe  an  amount  of  "weighting"  that  should  be  used  to  affect  a  particular  vertex  from 
the  skinCoord  field  of  the  Humanoid  node.  Each  item  in  this  list  has  a  corresponding 
index  value  in  the  skinCoordlndex  field  of  the  Joint  node,  which  indicates  exactly  which 
coordinate  is  to  be  influenced. 
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The  ulimit  and  llimit  fields  of  the  Joint  PROTO  speeify  the  upper  and 
lower  joint  rotation  limits.  Both  fields  are  three-element  MFFloats  eontaining  separate 
values  for  the  X,  Y  and  Z  rotation  limits.  The  ulimit  field  stores  the  upper  (i.e.  maximum) 
values  for  rotation  around  the  X,  Y  and  Z  axes.  The  llimit  field  stores  the  lower  (i.e. 
minimum)  values  for  rotation  around  those  axes.  Note  that  the  default  values  for  each  of 
these  fields  is  [],  which  means  that  the  joint  is  assumed  to  be  unconstrained. 

The  limitOrientation  exposedField  gives  the  orientation  of  the  coordinate 
frame  in  which  the  ulimit  and  llimit  values  are  to  be  interpreted.  The  limitOrientation 
describes  the  orientation  of  a  local  coordinate  frame,  relative  to  the  Joint  center  position 
described  by  the  center  exposedField. 

The  stiffness  exposedField,  if  present,  contains  values  ranging  between  0.0 
and  1.0  which  give  the  inverse  kinematics  system  hints  about  the  "willingness"  of  a  joint 
to  move  in  a  particular  degree  of  freedom.  For  example,  a  Joint  node's  stiffness  can  be 
used  in  an  arm  joint  chain  to  give  preference  to  moving  the  left  wrist  and  left  elbow  over 
moving  the  left  shoulder,  or  it  can  be  used  within  a  single  Joint  node  with  multiple 
degrees  of  freedom  to  give  preference  to  individual  degrees  of  freedom.  The  larger  the 
stiffness  value,  the  more  the  joint  will  resist  movement  in  the  corresponding  axis  (X,  Y, 
or  Z  for  entries  0,  1  and  2  of  the  stiffness  MFFloat  field). 

The  name  field  is  used  for  identifying  the  joints  at  runtime.  Each  Joint 
node  should  have  a  DEF  name  that  matches  the  name  field  for  that  Joint,  but  with  a 
distinguishing  prefix  in  front  of  it.  That  prefix  can  be  anything,  but  must  be  the  same  for 
all  the  Joint  nodes  in  a  particular  humanoid.  The  distinguishing  prefix  is  useful  in  the  case 
of  static  routing  to  the  Joint  nodes  of  multiple  humanoids  in  the  same  file.  If  only  a  single 
humanoid  is  stored  in  a  file,  the  prefix  should  be  "hanim"  (for  Humanoid  Animation). 
For  example,  the  left  shoulder  would  have  a  DEF  name  of  "hanim  l  shoulder".  The 
DEF  name  is  used  for  static  routing,  which  would  typically  connect 
Orientationinterpolators  in  the  humanoid  file  to  the  joints.  The  name  field  must  be 
present,  so  that  applications  are  able  to  identify  the  Joint  node  at  runtime. 

It  will  occasionally  be  useful  for  the  person  creating  a  humanoid  to  be  able 
to  add  additional  joints  to  the  body.  The  body  remains  humanoid  in  form,  and  is  still 
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generally  expeeted  to  have  the  basie  joints.  However,  they  may  be  thought  of  as  a 
minimal  set  to  whieh  extensions  may  be  added  (sueh  as  a  prehensile  tail).  See  the  seetion 
on  Non-standard  Joints  and  Segments.  If  neeessary,  some  of  the  joints  (sueh  as  the  many 
vertebrae)  may  be  omitted.  See  Appendix  A  for  details. 

3.  Deformation  Of  H-Anim  Body 

Realistie  representation  of  human  body  ean  only  be  aehieved  by  deformable 
humanoids.  H-Anim  2001  working  group  eonstrueted  a  mesh  based  deformation  proeess 
for  H-Anim  bodies.  Skin  tag  in  the  Humanoid  node  eontains  all  the  vertiees  those  build 
human  appearanee.  This  provides  integrity  of  the  skin  surfaee  between  the  joints.  By 
index  based  surfaee  building  proeess,  skin  eoordinates  may  be  read  from  skinCoord 
Index  tag.  Also  skin  eoordinate  normals  are  eontained  in  the  skinNormal  tag.  Eaeh  joint 
eontains  the  indiees  of  vertiees  in  skinCoordIndex  tag  .  This  list  of  indiees  will  be 
oriented  by  rotation  transformations  of  joint  they  belong  to  and  as  well  as  the  neighboring 
joints.  SkinCoordWeight  is  the  multiplieative  value  of  an  index.  This  value  shows  the 
pereentage  of  the  rotation,  whieh  will  be  applied  to  index 

This  design  provides  deformation  engines  all  information  they  need.  Following  is 
a  simple  deformation  engine  proeess  diagram. 
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Update  Rotation  value  of  connected  joint 


N 


Update  H-Anim  coordinates 

•  Read  joint 

•  Read  index 

•  Apply  rotation 

•  Update 


Figure  14.  Deformation  Process 


c. 


CONSTRUCTION  OF  VIRTUAL  HUMAN  FROM  LASER  SCAN 


1.  Method 

In  order  to  design  a  virtual  Human,  lasers  seans  from  Cyberware©  are  used. 
This  type  of  sean  gives  us  more  detailed  skin  information  than  the  other  seanning 
methods .  [DUTTON200 1  ] . 

a.  Data  Collection 

Colleetion  and  ordering  of  the  laser  sean  data  is  the  primary  task  in 
designing  a  virtual  human.  Every  seanning  system  has  its  own  file  format.  This  makes 
this  work  a  little  difficult.  However,  H-Anim  2001  standard  uses  a  skin  node,  which 
collects  all  vertex  information  and  constructs  a  surface  by  using  indices  from  this  field. 
By  direct  extraction  of  data  cloud  and  putting  it  into  this  field  provides  an  easy  solution. 
(Figure  15) 


51 


Figure  15.  Laser  Scan  Data  Cloud  Courtesy  of  Cyberware© 


Cyberware  Laser  scans  include  the  polygon  information  for  the  scanned 


body  skin.  This  information  is  used  for  mesh  construction.  This  body  mesh  has  a  great 
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number  of  triangles  in  order  to  provide  a  detailed  appearance. (Figure  16)  By  using  a 
triangle-reducing  program,  number  of  triangles  can  be  adjusted  for  performance  gain.  In 
this  research  a  free  version  of  3D  reducer©  from  Sergei  A  Zavalesky  is  used.  Triangle 
reduction  must  be  made  before  segmentation  process.  (Figure  17)  H-Anim  standard  does 
have  segment  node  fields  which  are  impossible  to  optimize  with  a  standard  tool. 


Figure  16.  Laser  Scan  Polygon  Mesh 


b.  Segmenting  Process 

For  segmentation  a  simple  bounded  cylinder  method  is  used.  Segmenter 
program,  written  in  Java  and  Virtual  Reality  Modeling  Language  (VRML),  performs  two 
tasks.  First  the  program  takes  segment  center  coordinates  for  each  body  part  as  user  input 
and  finds  all  possible  vertices  inside  a  cylinder.  Cylinder  coordinates  are  the  center  of 
two  segments,  which  are,  connected  each  other.  The  user  also  defines  end  points  of  end 
effectors.  The  second  task  is  writing  the  code  in  proper  VRML  format. 
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This  method  is  an  easy  solution  for  segmentation  proeess.  Previous  work 
[DUTTON2001]  used  a  eommereial  produet  for  this  purpose.  But  this  solution  is 
expensive  and  time  eonsuming.  Beeause  it  is  impossible  to  install  that  kind  of 
professional  modeling  programs  to  everyone’s  eomputer  and  moreover  modeling 
programs  have  a  long  learning  eurve.  But  modeling  programs  have  high  segmentation 
aeouracy  for  polygonal  mesh  and  more  eomplex  tools  for  rendering  proeess.  Aeeuraey 
near  end  points  and  around  joints  has  a  signifieant  effeet  for  deformation.  This  effeet  will 
be  diseussed  in  next  seetion. 

c.  Deformation  Issues 

Simple  Deformation  Engine  is  a  java  thread,  whieh  works  as  a  real-time 
rendering  engine  with  VRML  browser.  Its  performanee  is  direetly  proportional  to  the 
number  of  vertiees  inside  the  humanoid  mesh.  If  there  are  more  then  one  humanoids 
inside  the  virtual  environment  then  to  inerease  performanee  we  need  to  ehange  level  of 
detail  during  the  eonstruetion  of  humanoid  mesh. 

Another  issue  is  esthetie  appearanee  of  virtual  human.  Espeeially 
segmenting  through  false  joint  eenters  results  abnormal  artieulation  of  body  posture.  This 
ean  be  prevented  by  two  methods.  The  first  method  is  using  a  more  aeeurate  modeling 
tool.  Beeause  eaeh  person  will  have  his  own  humanoid  designed  onee.  So  someone  ean 
design  it  for  you  with  an  artistie  ability  and  aeouraey.  The  seeond  solution  is  defining 
joints  during  scanning  process.  Which  may  be  possible  by  putting  a  marker  to  the  joints. 
Another  way  is  standing  with  a  special  posture  that  shows  all  joint  tearing  points  during 
scanning. 

2.  Results 

Eigure  17  shows  the  final  product  of  the  segmentation  process.  Different 
articulations  are  experimented  with  different  orientation  around  x-  axis.  Eigures  18  and 
19  are  arm  and  leg  deformation  examples.  But  simple  deformation  engine  can  not  handle 
extreme  orientation  angles  as  seen  Eigure  20. 
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Figure  17.  Final  Product  of  Segmenter  Software 
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hanim_l_elbow 
x=0.917981 1924695969 
y=0 

z=0.{)2380006480962038 
angle=-0.5262548 1 03928566 

X  y  z  angle 


Figure  18.  Arm  Deformation  Example 
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hanimiknee 

x=0.6928366608917713 

y=0 

z=-0.027364 1 1988735199 
angle=0.66342959 1 8047428 


X  y  z  angle 


Figure  19.  Leg  Deformation  Example 
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Figure  20.  Deformation  Problems 

Human  skin  deformation  is  very  eomplex  process.  Muscles  under  the  skin  have 
nonlinear  stretching  and  expanding  structure.  Even  though  the  joint  node  contains  a 
weight  field,  which  provides  deformation  coefficient  for  orientation,  writing  a  complex 
deformation  engine  may  decrease  the  performance.  That’s  why  we  have  static  values  for 
this  field.  And  static  values  may  create  abnormal  stretching  polygons  over  the  surfaces 
around  the  joints  as  seen  in  Figure  20. 
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IV.  VIRTUAL  HUMAN  INTERFACE  FOR  MARG  BODY 

TRACKING  SYSTEM 


A.  ANALYSIS  OF  MARG  HARDWARE 

1,  Orientation  Signals 

As  stated  before  MARG  system  uses  three-axis  information  from  three  sensors. 
Magnetic  and  accelerometer  sensor  give  correction  for  drift  and  rate  sensor  give 
orientation  information  after  using  correction.  All  these  nine  signals  are  reduced  to  four 
tuples  which  are  unit  quaternion  and  defined  as  rotation  around  given  vector. 
[BACHMANN2000]  explains  reduction  of  the  X-matrix  in  detail.  Optimized  virtual 
human  has  fifteen  segments  as  shown  in  Figure  21. 


Figure  2 1 .  Segments  for  a  virtual  human  in  MARG 

Currently  there  does  not  exist  any  graphics  rendering  hardware,  which  uses 
directly  quaternion  information.  That  is  the  reason  why  this  system  uses  VRML 
language.  VRML  language  rotations  are  defined  as  axis  rotations,  which  is  very  similar  to 
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quaternion  rotations.  That’s  makes  it  easy  to  convert  quaternion  rotations  to  axis  angle 
rotations.  The  drawback  of  this  conversion  is  singularity  of  some  angles. 

New  programmable  vector  graphics  engines  may  provide  a  quaternion  based 
calculation  system.  This  approach  will  increase  the  performance  of  the  system. 

Our  system  uses  quaternion  orientation  information  over  network.  So  network 
speed  will  directly  affect  frame  rate.  New  algorithms,  which  may  reckon  articulations 
during  loss  of  network  information,  may  be  developed  as  a  new  research  project. 

2.  Location  Signals 

To  locate  a  human  body  this  research  used  GPS  signals.  Today  there  exist  very 
small  sized  GPS  instruments,  which  can  directly  be  plugged  to  wearable  computers.  A 
parser  class  is  written  to  read  this  information  from  a  GPS  by  using  NMEA  standard 
183.  Location  signals  are  taken  as  Latitude  and  Longitude  and  applied  to  the  human  root 
translation  field.  Because  this  field  is  the  root  of  all  joints  and  can  translate  all  segments 
at  the  same  time. 

Error  rates  applied  to  GPS  signals  by  government  sometimes  make  it  impossible 
to  use  them  for  human  body  location.  That’s  why  we  used  a  differentiation  method  for 
localization.  Lor  this  purpose  two  GPS  instruments,  one  for  tracker  one  for  the  actual 
person  whose  body  is  being  tracked,  are  used.  The  difference  between  these  two  locations 
is  taken  and  the  location  of  the  tracker  is  used  as  a  reference  to  locate  the  actual  tracked 
person.  The  error  for  the  tracker  and  tracked  person  is  same  for  the  close  regions.  So 
relative  position  will  be  as  accurate  as  actual  location.  (Ligure  22)  In  other  words  PtPr 
will  parallel  to  PtgpsPrgps-  For  the  implementation  this  vector  will  be  brought  onto  a 
special  origin  defined  by  the  user. 
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Figure  22.  GPS  Positioning  by  two  GPS  Instruments 

B.  IMPLEMENTATION  AND  RESULTS 

Figure  23  is  general  architecture  of  the  implemented  system.  Conceptually  the 
system  is  designed  at  least  for  fifteen  MARG  sensors.  But  currently  there  exists  only  one 
sensor  for  experiments.  However  designing  one  channel  is  enough  for  the  proof  of 
concept  because  all  other  fourteen  channels  will  have  same  architecture. 
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F igure  2  3 .  General  Architecture 

MARG  Research  group  is  working  on  wireless  system,  which  is  not  implemented 
yet.  Another  part  of  the  project  deals  with  networking  process.  Because  of  these  two 
ongoing  researches  this  work  implemented  a  simulated  network  server  to  read  the 
quaternion  orientation  information.  Simulated  server  sends  previously  recorded  data 
through  network  to  client  software.  Motion  capture  is  done  by  MARG  system  designed 
by  [BACHMANN2000]. 

Other  function  of  Server  is  sending  location  information.  Commercial  GPS 
connected  to  serial  port  of  computer  updates  position  field  of  server.  Special  parser  thread 
written  in  JAVA  language  reads  communication  port  of  computer,  and  decodes  NMEA 
183  standard  data  into  text  format.  Finally  the  server  program  packs  this  data  with 
quaternions. 
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There  are  two  possible  package  types  for  this  system.  One  of  them  is  packing  all 
rotation  info  into  one  package.  Other  possibility  is  sending  this  rotation  info  separately 
for  each  joint.  Implemented  system  used  second  packet  format  as  shown  in  Figure  24. 


UDP  HEADER 


SEGMENT  NAME 


X  YZ  W 


Eigure  24.  Packet  Eormat 

Client  program  reads  these  network  packages  and  updates  the  joint  rotations  and 
the  humanoid  root  translation  field.  During  this  process  simple  deformation  thread 
updates  the  skin  in  order  to  create  seamless  body  form. 

Performance  of  the  system  as  frame  per  second  is  shown  in  Eigure  24.  A  Pentium 
IV  with  GEORGE  2.5,  32  Mb  graphics  card  has  been  used  to  make  the  experiments. 
Clearly  level  of  detail  is  main  reason  for  the  decrease  in  performance. 
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Figure  25.  Performance  Analysis 


Another  parameter  effecting  to  the  system  performance  is  network  delays. 
Network  latency  represents  one  of  the  biggest  challenges  for  net-VE  program 
writers.  [ZYDA2000].  In  order  to  decrease  this  effect  server  program  sends  user  datagram 
packets  (UDP)  to  the  network.  Disadvantage  of  the  UDP  is  that  they  are  not  reliable.  But 
number  of  packets  sent  for  a  human  body  is  so  many  that  loss  of  some  packets  will  not 
directly  affect  the  appearance. 

Simple  deformation  engine  is  the  most  important  and  time-consuming  thread  for 
the  rendering  performance.  Depending  of  the  number  of  vertices  this  engine  renders  the 
surface  data  for  each  time  a  new  rotation  quaternion  package  updates  the  humanoid  joint 
structure. 
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V.  SUMMARY  AND  CONCLUSION 


A.  SUMMARY 

This  thesis  presented  a  real-time  method  for  the  speeifieation  of  seamless  virtual 
human  interfaee  for  MARG  body  traeking  system  based  on  H-Anim  2.0  standard. 
Implementation  is  realized  by  Java- VRML  language. 

Designing  a  virtual  human  has  three  aspeets.  First  aspeet  is  eonstruetion  of  a 
realistie  graphieal  representation  of  human.  This  thesis  used  Cyberware©  laser  seans  in 
order  to  aehieve  this  goal.  Seeond  aspeet  is  related  to  animation  of  this  statie  strueture. 
For  this  purposes,  laser  sean  data  eloud  is  segmented  into  a  skeleton  strueture.  This 
skeleton  strueture  is  a  low  level  of  detailed  version  of  H-Anim  2.0  speeifieation.  H- 
Anim  2.0  standard  also  provided  a  seamless  human  skin  strueture. 

Final  aspeet  of  virtual  human  design  is  to  ereate  a  realistie  and  real-time 
artieulation  for  eaeh  segment.  Soureeless  MARG  sensor  system,  whieh  is  designed  by 
Naval  Postgraduate  Sehool  researehers,  is  suitable  motions  eapture  system  for  this 
purpose.  It  uses  singularity-free  quaternions  for  rotation.  Quaternions  provide  not  only  a 
safe  mathematieal  base  for  rotations  but  also  speed  for  the  ealeulations  of  orientations. 

Motion  eapture  usually  deals  only  with  segment  rotations.  This  researeh  added 
position  traeking  to  aehieve  the  goal  for  a  remote  body  traeking.  Global  positioning 
system  is  regarded  as  an  optimal,  eheap  and  global  method  to  loeate  the  human  body. 


B,  CONCLUSIONS 

Creating  a  realistie  virtual  world  ean  be  measured  primarily  by  human  pereeption. 
The  more  natural  the  objeets  are,  the  more  eomfortable  beeomes  the  human  in  this 
environment.  But  this  also  brings  a  problem.  Humans  beeome  more  seleetive  in  realistie 
virtual  environment.  As  their  expeetation  for  reality  inereases,  their  forgiveness  for  any 
type  of  error  deereases. 

Keeping  this  in  mind  this  researeh  has  following  eonelusions. 
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Realistic  human  body  can  be  created  by  directly  using  laser  scans  and  segmenting 
this  scan  by  either  a  commercial  3-D  modeling  software  or  a  simple  segmentation  VRML 
program. 

Segmented  human  model  will  be  able  to  get  information  for  each  segment  through 
network  by  using  MARG  body  tracking  system. 

Location  information  is  usually  regarded  as  secondary  issue  for  body  tracking. 
This  research  takes  this  issue  as  one  of  primary  problems  and  implemented  a  solution  by 
using  GPS  parser.  But  the  error  in  GPS  signals  and  lack  of  navigational  position 
information  in  virtual  worlds  create  drawback.  Differential  GPS  systems  can  prevent  this 
disadvantage  by  using  different  techniques. 

Speed  of  Internet  and  increase  of  CPU  power  enabled  ordinary  users  to  be  a  part 
of  3D  virtual  worlds.  Human  existence  in  these  worlds  is  currently  possible  by  static 
avatars.  Some  of  these  avatars  have  limited  animated  behaviors  for  simple  body 
movements  like  hand  waving.  Real-time  sourceless  body  tracking  will  make  it  possible 
for  these  avatars  to  have  infinitely  many  articulations  with  realistic  skin  deformation. 
This  research  proved  this  possibility  by  using  MARG  body  tracking  system  which  is, 
applied  a  VRML-H-Anim  based  avatar  with  GPS  support. 


66 


LIST  OF  REFERENCES 


BACHMANN99 


BACHMANN2000 


CRAI89 


COOKE92 


DUMAN99 


DURL95 


DUTTON2001 


Bachmann,  E.,  Duman  L,  EFsta,  EF.,  McGhee  R.,  Yun,  X.,  and  Zyda, 
M.,  “Orientation  traeking  for  Humans  and  Robots  Elsing  Inertial 
Sensors,”  International  Symposium  on  Computational  Intelligence 
in  Robotics  <&  Automation  (CIRA  99),  Monterey,  CA,  November 
1999,  pp.187-194. 

Bachmann  Eric  ,  Inertial  and  MagneticTtracking  of  Limb  Segment 
Orientation  Inserting  Humans  Into  Synthetic  Environments, 
Dissertation  ,  Naval  Postgraduate  School,  Monterey,  CA,  March 
2000 

Craig,  J.,  Introduction  to  Robotics:  Mechanics  and  Control,  Second 
Edition,  Addison- Wesley  Publishing  company,  Inc.,  Menlo  Park, 
CA,  1989. 

Cooke,  J.,  Zyda,  M.,  Pratt,  D.,  McGhee,  R.,  “NPSNET;  Plight 
Simulation  Modeling  Using  Quaternions,”  Presence,  Vol.  1,  No.  4, 
MIT  Press,  Cambridge  MA,  Pall,  1992,  pp.  404-420. 

Duman,  I.,  Design,  Implementation  and  Testing  of  a  Real-Time 
Software  System  for  a  Quaternion-Based  Attitude  Estimation  Eilter, 
Master’s  Thesis,  Naval  Postgraduate  Sehool,  Monterey,  CA,  March 
1999. 

Durlach,  N.,  and  Mavor,  A.,  National  Research  Council,  Virtual 
Reality:  Scientific  and  Technological  Challenges,  National 
Academy  Press,  Washington,  DC,  1995. 

Dutton  James  Allen  ,  Developing  Articulated  Human  Models  Erom 
Laser  Scan  Data  For  Use  As  Avatars  In  Real-time  Networked 
Virtual  Environments  ,  Master’s  Thesis,  Naval  Postgraduate 
School,  Monterey,  CA,  March  2001 


67 


EBERLY99 

HANIM2001 

FRAN69 

EOEE97 

HAND93 

HAMILTON  1899 

HODGINS2000 

MCGHEE93 

MCGHEE97 

MORSE53 


David  Eberly  Quaternion  Algebra  and  Ca/cM/w5,  Magic 
Software  6006  Meadow  Run  Court  Chapel  Hill,  NC  27516, 
eberlv@maqic-software.com 

H-Anim  organization,  H-Anim  2.0  draft  specification,  2001 
www.hanim.org 

Frank,  A.  A.,  McGhee,  R.  B.,  Some  Considerations  Relating  to  the 
Design  of  Autopilots  for  Legged  Vehicles,  Journal  of 
Terramechanics,  Vol.  6,  No.l,  pp.  23-35,  Pergamon  Press,  Great 
Britain,  1969. 

Foley,  J.  D.,  Van  Dam,  A.,  Feiner,  S.  K.  and  Hughes  J.  F., 
Computer  Graphics  Principles  and  Practice,  Second  Edition  in  C, 
Addison- Wesley,  1997. 

Hand,  C.,  “A  Survey  of  3-D  Input  Deviees'”,  Technical  Report  CS 
TR94/2,  De  Montfort  University,  Leicester,  UK,  1993. 

Hamilton  Sir  William  R.  ,  On  quaternions  ,Royal  Irish  Academy 
press  1899 

Hodgins  Jesica,  Animating  Human  motion  Seientifie  American 
1998 

http://www.sciam.com/1998/0398issue/0398hodgins.html 

MeGhee,  R.  B.,  CS4314  Class  Notes:  Derivation  of  Body  Angular 
Rates  to  Euler  Angle  Rates  Relationship,  Naval  Postgraduate 
Sehool,  Monterey,  CA,  1993. 

MeGhee,  R.  B.,  CS4314  Class  Notes:  Finite  and  Infinitesimal 
Quaternion  Rotations,  Naval  Postgraduate  Sehool,  Monterey,  CA, 
1997. 

Morse,  Philip  M.;  Feshbach,  Herman  Methods  of  Theoretical 
Physics  McGraw-Hill  1953  . 


68 


SHOEMAKE85 

SHOEMAKE87 

SKOP96 

THAEMANN97 

USTA99 

WATT98 

WIEKINS99 

YUN99 

ZYDA99 


Shoemake,  K.,  “Animating  Rotation  with  Quaternion  Curves”, 
SIGGRAPH85,  pp.  245-254,  ACM  Press,  1985. 

Ken  Shoemake,  Animating  rotation  with  quaternion  calculus, 
ACM  SIGGRAPH  1987,  Course  Notes  10,  Computer  Animation; 
3-D  Motion,  Specifieation,  and  Control. 

Skopowski,  P.,  Immersive  Articulation  of  the  Human  Upper  Body 
in  a  Virtual  Environment,  Master’s  Thesis,  Naval  Postgraduate 
Sehool,  Monterey,  CA,  Deeember  1996. 

Daniel  Thalmann  ,  Virtual  Humans  in  Virtual  Environments  A  New 
View  of  Multimedia  Applications  1997  Computer  Graphics  Eab 
(EIG),  Swiss  Eederal  Institute  of  Technology  (EPEE)  CH-1015 
Eausanne,  Switzerland 

Usta,  U.,  Comparison  of  Quaternion  and  Euler  Angle  Methods  for 
Joint  Angle  Animation  of  Human  Figure  Models,  Master’s  Thesis, 
Naval  Postgraduate  School,Monterey,  CA,  March  1999. 

Watt,  A.,  and  Watt,  M.,  Advanced  Animation  and  Rendering 
Techniques,  Theory  and  PracticQ,  ACM  Press,  New  York,  NY, 
1992. 

Wilkins  ,  David  R.  On  Quaternions  from  Sir  William  Hamilton 
Trinity  College,  Dublin  1999 

http://www.maths.tcd.ie/pub/HistMath/People/Hamilton/ 

Yun,  X.,  Bachmann,  E.,  McGhee  R.,  Whalen,  R.,  Roberts,  R., 
Knapp,  R.,  Healey,  A.,  and  Zyda  M., “Testing  and  Evaluation  of  an 
Integrated  GPS/INS  System  for  Small  AUV  Navigation,”  IEEE 
Journal  of  Ocean  Engineering,  Vol.  24,  No.  3,  July  1999.  pp.  396- 
404. 

Zyda  M.,  and  Singhal,  S.,  Networked  Virtual  Environments,  Design 
and  Implementation,  ACM  Press,  New  York,  NY,  1999. 


69 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


70 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 


2.  Dudley  Knox  Library 

Naval  Postgraduate  School 
Monterey,  California 


3.  Alper  SINAV 
Turkish  Navy 
Istanbul  Turkey 

4.  Prof. Michael  Zyda 
Naval  Postgraduate  School 
Monterey  CA  93943-5101 


5.  Prof.  Beny  Neta 

Naval  Postgraduate  School 
Monterey  CA  93940-5101 


6.  Prof.  Xiaoping  Yun 
Naval  Postgraduate  School 
Monterey  CA  93940-5101 

7.  Deniz  Kuvvetleri  Komutanligi 
Ankara  Turkiye 


71 


